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CHAPTER  I 


INTRODUCTION 

This  dissertation  develops  a  technique  for  optimizing  a 
spares  provisioning  problem  over  a  multiyear  period.  The  spares 
provisioning  problem  has  been  previously  formulated  [1]  and  a 
heuristic  optimization  procedure  specified.  The  complexity  of  the 
problem  is  such  that  it  could  not  be  solved  exactly  by  standard 
techniques  at  that  time.  Through  the  use  of  an  amalgamation  of 
branch  and  bound,  partial  enumeration,  and  dynamic  programming 
techniques,  an  approach  that  assures  convergence  to  the  optimal 
solution  is  obtained  here.  This  work  discusses  and  justifies  the 
solution  technique  to  the  spares  provisioning  problem. 

1.1  Background 

The  spares  provisioning  problem  is  a  stochastic  problem  in 
which  an  optimum  policy  for  purchasing  spare  equipment  and  repair 
channels  over  a  multitime  period  horizon  is  sought.  For  each  time 
period,  a  specified  number  of  machines  are  required  to  be  on  line. 

As  normally  happens,  some  of  these  machines  will  break  down.  Upon 
breakdown,  a  machine  is  removed  to  a  repair  station  and  replaced  by 
a  spare,  if  one  is  available.  To  maintain  a  reasonable  level  of 
operation  there  is  a  requirement  placed  upon  the  system  that  a  spare 
be  available  a  specified  percentage  of  the  time. 

There  is  a  tradeoff  to  be  made  between  the  number  of  spares 
purchased  and  the  number  of  repair  facilities  set  up.  If  a  machine 
must  wait  to  be  repaired  because  all  of  the  repair  facilities  arc 
busy,  that  machine  is  not  available  for  use  as  a  spare  for  a  longer 
period  of  time.  The  greater  the  number  of  repair  facilities,  the 


greater  the  likelihood  that  a  machine  will  not  wait  to  be  repaired. 

On  the  other  hand,  if  there  is  an  abundance  of  spare  machines  avail¬ 
able,  a  delay  in  the  repair  of  the  machine  may  not  impact  on  the 
system.  It  is  desired,  therefore,  depending  on  the  cost  of  the 
spares,  the  cost  of  setting  up  a  repair  facility,  and  the  relative 
rates  of  machine  failure  and  repair,  to  determine  the  most  cost 
effective  mix  of  spares  and  repair  facilities  that  satisfies  the 
availability  requirement. 

The  problem  is  further  complicated  because  the  optimization 
is  over  multiple  time  periods.  At  the  beginning  of  each  time  period 
the  number  of  machines  on  line  may  be  changed.  Either  new  machines 
may  be  added  or  old  machines  removed.  Because  this  problem  is  over 
a  multiperiod  horizon,  the  optimal  policy  of  selecting  spares  or 
repair  channels  for  each  individual  time  period  is  not  necessarily 
the  optimal  policy  when  considered  in  the  context  of  the  total  hori¬ 
zon.  Repair  channels  and  spares  that  may  not  be  needed  in  the 
current  time  period  but  that  will  be  required  in  a  future  time  period 
may  be  purchased  in  the  current  period  because  of  lower  cost. 

From  the  above  description  of  the  problem  it  might  appear 
that  obtaining  a  solution  should  not  be  too  difficult.  It  is  a 
multiperiod  problem  for  which  there  are  two  independent  variables 
(number  of  spares  and  number  of  repair  channels)  for  each  time 
period.  There  is  also  a  single  constraint  imposed  on  the  avail¬ 
ability  of  replacement  machines  for  each  of  the  time  periods  and  an 
objective  function  associated  with  the  purchase  cost  of  spares  and 
machines.  If  the  objective  and  constraining  functions  had  reason¬ 
ably  simple  structures,  the  problem  could  be  solved  using  standard 
integer  linear  or  nonlinear  solution  techniques.  The  difficulty, 
however,  is  that  the  constraint  equations  do  not  have  a  nice  func¬ 
tional  dependence  on  the  independent  variables.  Because  of  the 
stochastic  nature  of  the  problem,  there  is  a  complicated  dependence 
on  the  variables. 

The  underlying  stochastic  process  can  be  formulated  in  terms 
of  the  classical  machine  repair  problem.  The  work  of  Taylor  and 
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Jackson  [2]  was  one  of  the  earliest  applications  of  queueing  theory 
to  the  spares  provisioning  problem.  In  this  formulation,  the  ser¬ 
vice  time  and  the  time  between  failures  are  assumed  to  be  exponen¬ 
tially  distributed  random  variables.  As  will  be  shown  subsequently, 
the  steady  state  probability  that  n  machines  are  in  need  of  repair 
can  be  specified  as  a  function  of  the  number  of  machines  in  the 
system,  the  number  of  repair  channels,  the  number  of  spares,  and 
the  relative  failure  rate  and  repair  rate.  The  explicit  formulation 
is  dependent  on  whether  the  number  of  repair  channels  is  greater  or 
less  than  the  number  of  spares  and  whether  n  ,  the  number  in 
repair,  is  greater  or  less  than  the  number  of  channels  or  number  of 
spares.  As  will  be  seen,  it  is  this  aspect  of  the  problem  that 
causes  the  greatest  complication  since  the  explicit  forms  of  the 
constraints  cannot  be  specified  until  the  policy  is  known. 

Another  facet  of  the  problem,  which  is  a  result  of  the  multi¬ 
ple  time  periods,  is  the  determination  of  the  machine  failure  rate. 

As  part  of  the  formulation  of  the  problem  (see  Reference  [1]),  the 
failure  rate  specified  for  each  time  period  is  the  average  failure 
rate  taken  over  the  number  of  machines  introduced  in  the  time  period, 
the  number  of  machines  repaired  in  the  previous  time  period,  and  the 
machines  from  the  previous  time  period  that  had  not  been  repaired. 

In  that  the  expected  number  of  machines  repaired  in  a  time  period 
is  dependent  on  the  number  of  repair  channels,  the  average  failure 
rate  for  a  given  time  period  is  dependent  on  the  number  of  repair 
channels  in  previous  time  periods. 

Therefore,  what  may  at  first  appear  to  be  a  simple  problem 
has  become  very  complex  in  that  the  functional  dependence  is  not 
explicitly  known  until  a  policy  is  enumerated,  and  a  key  parameter, 
the  failure  rate,  is  dependent  on  the  selected  policy  of  previous 
time  periods.  Because  of  these  two  complicating  factors  a  solution 
technique  that  combined  partial  enumeration  and  the  sequential 
technique  of  branch  and  bound  was  needed  to  solve  the  problem. 


4 


1.2  Purpose 

Spares  provisioning  analysis  has  been  applied  to  the 
maintenance  of  marine  gas  turbine  engines  for  the  U.S.  Department 
of  the  Navy  [1].  For  this  particular  problem  the  reliability  of 
a  new  gas  turbine  engine  is  assumed  to  depend  on  its  time  of  intro¬ 
duction.  The  reliability  of  a  used  item  is  assumed  to  depend  on 
the  time  of  its  last  repair.  It  is  assumed  that  there  is  a  learn¬ 
ing  curve,  and  that  with  increased  experience,  the  gas  turbine 
engines  will  have  an  improved  reliability.  Therefore,  the  failure 
rate  used  in  the  formulation  of  this  problem  is  dependent  on  the 
number  of  engines  introduced  in  each  time  period  as  well  as  on  the 
number  repaired  in  previous  time  periods.  Through  this  failure 
rate  the  availability  constraint  for  a  given  time  period  is  coupled 
to  the  policy  used  in  all  of  the  previous  time  periods.  This  kind 
of  situation  may  be  typical  of  many  repairable  item  inventory  sys¬ 
tems  in  which  the  items  to  be  controlled  result  from  new  tech¬ 
nology  . 

The  determination  of  an  optimal  purchase  policy  for  the  mix 
of  spares  and  repair  channels  needed  in  each  time  period  was  pre¬ 
viously  attempted  using  a  heuristic  approach  which  separated  the 
problem  by  time  periods  [1].  While  this  approach  provided  an 
apparently  "good"  policy — and  a  feasible  one  in  terms  of  its  satis¬ 
faction  of  the  spares  availability  constraint — there  was  no  assur¬ 
ance  that  it  was  optimal,  and  it  has  been  shown  to  be  sub-optimal 
for  some  cases. 

This  dissertation  employs  the  same  basic  model  as  was  used 
in  [1],  but  approaches  the  optimization  in  a  different  manner. 

The  reliability  for  each  time  period  is  initially  approximated  by 
ignoring  the  terms  dependent  on  the  previously  repaired  components, 
and  an  optimum  policy  for  the  approximate  problem  is  determined  over 
all  time  periods.  The  resulting  trial  solution  is  checked  for 
feasibility.  If  the  solution  is  not  feasible,  branching  occurs. 
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creating  a  new  set  of  reliability  factors.  The  new  set  of  prob¬ 
lems,  using  the  updated  reliability  factors,  more  closely  bounds 
the  original  problem.  This  branch  and  bound  technique  proceeds 
until  the  optimal  solution  is  found. 

The  use  of  this  approach  still  results  in  an  integer  pro¬ 
gramming  model.  Through  techniques  discussed  within  this  paper, 
the  lattice  of  integers  that  define  the  boundaries  of  the  constraint 
set  are  determined.  By  taking  the  convex  hull  of  this  constraint 
set  and  using  a  linear  objective  function,  the  problem  could  then 
be  solved  using  the  simplex  method.  However,  it  was  found  to  be 
computationally  simpler  to  solve  the  problem  using  dynamic  program¬ 
ming. 


The  technique  described  in  this  dissertation  not  only  gives 
an  optimal  solution  to  the  spares  provisioning  problem,  but  also 
requires  less  computer  capability  and  is  computationally  faster 
than  the  previous  heuristic  approach.  Furthermore,  the  method  is 
applicable  to  a  more  general  type  of  problem  and  can  be  used  to 
solve  the  particular  spares  provisioning  problem  with  a  more 
complex  cost  function. 


1.3  Overview  of  the  Study 

This  work  consists  of  four  additional  chapters  plus  two 
appendices.  Chapter  II  presents  the  spares  provisioning  problem 
in  detail,  with  an  explicit  mathematical  characterization.  The 
heuristic  solution  suggested  in  [1]  is  also  discussed.  In  Chapter 
III  the  mathematical  development  of  the  optimization  techniques 
is  discussed.  This  includes  a  discussion  of  the  convergence  of  the 
solution,  the  technique  for  the  generation  of  the  constraint  set, 
and  the  dynamic  programming  solution.  Computational  results  and 
example  problems  with  solutions  are  presented  in  Chapter  IV.  The 
final  chapter  summarizes  this  study  and  presents  conclusions  and 
recommendations . 


CHAPTER  II 


PROBLEM  DEFINITION 

While  most  of  the  effort  described  in  this  dissertation  has 
dealt  with  the  solution  of  a  particular  problem,  an  optimal  spares 
provisioning  policy  for  maintaining  an  adequate  number  of  operating 
marine  gas  turbine  engines,  the  techniques  developed  herein  can  be 
applied  to  a  more  general  class  of  problems.  This  chapter  discusses 
the  particular  gas  turbine  engine  problem,  its  formulation,  and  a 
previously  specified  heuristic  approach  for  obtaining  a  good  policy. 

II. 1  The  Gas  Turbine  Engine  Problem 

The  immediate  problem  addressed  here  is  one  of  determining 
the  minimum  cost  policy  for  purchasing  spares  and  repair  channels 
that  assures  that  there  is  an  adequate  supply  of  spares  available 
for  replacing  marine  gas  turbine  engines*  that  randomly  fail.  A 
multiyear  planning  horizon  is  considered,  with  the  possibility  that 
additional  turbines  may  be  introduced  in  each  year,  requiring  addi¬ 
tional  spares  and/or  repair  channels. 

When  a  turbine  fails,  a  request  is  made  for  a  replacement. 
This  replacement  comes  from  the  pool  of  spares,  and  the  component 
that  failed  is  removed  and  taken  to  the  repair  facilities.  If  at 
least  one  repair  channel  is  not  busy,  repair  begins  immediately; 


*There  are  two  components  in  a  gas  turbine:  the  gas  gener¬ 
ator  and  the  power  turbine.  Each  component  type  requires  its  own 
individual  set  of  spares  and  repair  channels,  and  thus  we  can  con¬ 
sider  independent  provisioning  systems.  We  will  use  the  term  "tur¬ 
bine,"  but  in  reality  we  are  dealing  with  either  the  gas  generator 
or  the  power  turbine  components. 
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» 

if  all  facilities  are  busy,  the  turbine  joins  the  queue  and  waits 
for  service.  Upon  completion  of  the  repairs,  the  gas  turbine  is 
placed  into  the  spares  pool.  If  there  is  an  inadequate  number  of 
repair  facilities,  a  queue  of  nonfunctioning  turbines  will  develop; 
as  the  queue  builds  and  the  spares  pool  diminishes,  it  is  possible 
for  the  spares  pool  to  become  depleted.  If  a  failure  occurs  when 
the  spares  pool  is  empty,  there  wil,l  be  no  available  replacement 
and  the  ship  in  which  the  failed  turbine  resides  will  not  be  able 
to  function.  Therefore,  there  is  imposed  on  the  problem  the  con¬ 
straint  that  the  probability  that  a  spare  is  available  when  a 
turbine  fails  be  at  least  0.90  (or  any  other  value  specified). 

Through  either  a  learning  process  or  a  conscious  component 
improvement  program  (CIP) ,  it  is  assumed  that  the  reliability  of  a 
turbine  improves  both  with  successive  years  of  manufacture,  and 
after  each  repair.  In  conjunction  with  these  assumptions  is  the 
assumption  that  the  mean  failure  rate  for  a  time  period  is  the 
average  of  the  failure  rate  for  the  new  machines  introduced  in  the 
time  period,  the  failure  rate  for  the  turbines  repaired  in  the  pre¬ 
vious  time  period,  and  the  mean  failure  rate  assumed  for  the  previ¬ 
ous  time  period  for  those  turbines  that  were  not  repaired  in  that 
time  period. 

This  problem  was  formulated  in  [1]  as  the  "classical" 
machine  repair  with  spares  queueing  problem.  Failure  and  service 
times  are  assumed  to  be  exponentially  distributed  random  variables. 
The  formulation  of  Reference  [lj  is  presented  here.  Based  upon  this 
framework,  a  mathematical  programming  problem  is  defined.  The  solu¬ 
tion  technique  for  this  problem  is  given  in  the  next  chapter. 

II. 2  Queueing  Model 

In  the  classical  machine  repair  with  spares  problem,  the 
machines  may  be  in  one  of  four  distinct  states.  They  may  be  up  and 
operating;  down  and  in  the  queue  awaiting  repair;  undergoing  repair; 
or  in  the  spares  pool,  where  they  are  available  as  replacements. 


« 
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For  a  given  set  of  parameters  (the  service  rate,  repair  rate,  number 
of  machines,  number  of  spares,  and  number  of  repair  channels)  under 
the  assumption  of  exponential  service  and  failure  distributions, 
the  system  will  eventually  reach  a  steady  state;  that  is,  the  prob¬ 
ability  that  a  machine  is  in  a  particular  state  is  independent  of 
time.  Since  the  steady  state  solution  of  the  classical  machine  re¬ 
pair  with  spares  problem  is  known  and  the  transient  model  appears 
to  be  analytically  intractable,  steady  state  is  assumed  in  this 
dissertation.  In  reality,  the  problem  never  reaches  steady  state 
since  new  gas  turbines  are  introduced  at  random  intervals  of  time 
throughout  the  year.  However,  it  is  assumed  that  new  turbines  are 
introduced  only  at  the  beginning  of  each  year  and  that  the  transient 
effects  caused  by  the  introduction  of  turbines  damp  out  immediately. 
Therefore,  each  time  period  analyzed  uses  the  steady  state  formula¬ 
tion. 


The  transient  effects  are  currently  under  investigation,  but 
since  this  dissertation  is  concerned  with  the  optimization  of  the 
purchase  policy,  the  steady  state  assumption  is  accepted  without 
further  discussion  [3]. 


To  understand  the  formulation  of  this  problem,  the  following 
notation  is  introduced. 


P 


n,  i 


steady-state  probability  that  n  items  are  "down" 
in  year  i 


c^  =  number  of  repair  facilities  to  plan  for  year  i 

y^  =  number  of  spares  to  plan  for  year  i 

=  component  failure  rate  (Poisson  mean)  for  new  or 
repaired  items  in  year  i 

Aj,  =  average  failure  rate  for  the  population  in  year  i 

R.  =  expected  number  of  components  repaired  for  year  i 

=  component  population  size  in  year  i 

1 /Mj  =  average  service  (removal,  transportation,  and  repair) 
time,  exponential  mean,  for  year  i  . 
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For  steady  state,  the  following  equations  are  obtained  [1] : 


n!  V^i /  0,1 


K  (\X 

n-c.  1  p.  )  p0,i 
c  1  c  !  V  V 

I  i  C±. 


(0<n<c^) 


(c^n^y^  b 


Pn4 


yi 

Mi  V 


(M^-n+y^) !  ci  ci 


/M" 


7 — -\t)  po,i  <vii"S'i+Mi> 


(II-l) 


n!  \yi/  P°«i 


(0<n<y . ) 


Mi\!  (\X 


(M^-n+y^) !  n!  I  ]rJ  p0,i  ^i^i*  <ci=yi>  e 


M^M.!  /  A.\n 

— “ — - 1  —  1  p„  .  (  c .  <n<y .  +M . ) 

n-c  \Vil  0»i  i-=-  i  i 


(M.-n+y.)!  c.  *  c.!  '  ^ 
iii  l  i 


M.+y  . 

P0,i  :  £0  Pn*i 


f(Mi-Mi-l)Xi  +  Ri-lXi-l  +  (Mi-l-Vl>Xi-li/Mi  *  Mi^Mi-l 


(f(Mi-Mi-l)Xi  +  Ri-lXi-l  +  <Mi-l-Ri 
(fRi-lXi-l+  (Mi-l  "  Ri-l)Xi-l}/Mi-l  ’ 


Mi^i 

?i  =  \  Mi  ■  L  (n-y^p  i 

1  n=y^-KL  1  n’x 


Mi^Mi-l 


(II-2) 


(II-3) 
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A  review  of  these  equations  gives  insight  into  the  nature  of 

the  problem.  Equation  (II-l)  gives  the  probability  that  the  system 

is  in  state  n  for  the  year  i  ;  a  state  is  defined  by  the  number 

of  machines  that  are  down.  The  form  of  the  expression  for  p  .  is 

n,  i 

dependent  on  the  relative  values  of  n^  ,  ,  and  y^  .  Equa¬ 

tion  (II-l)  consists  of  six  different  expressions.  The  first  three 
( a,b,a )  hold  if  there  are  at  least  as  many  (if  not  more)  spares  as 
there  are  repair  channels.  The  second  three  expressions  ( d,e,f) 
hold  if  the  number  of  repair  channels  is  greater  than  the  number 

of  spares.  The  exact  expression  for  p  .  is  dependent  on  whether 

n,i 

n  is  greater  than,  less  than,  or  equal  to  c^  or  y^  .  It  is 

apparent,  therefore,  that  it  is  not  possible  to  compute  the  state 
of  the  system,  as  determined  by  the  probability  that  n  machines 
are  down,  until  a  history  of  the  number  of  spares  and  repair  chan¬ 
nels  is  specified.  This  aspect  of  the  "undetermined  form  of 

p  .  "  is  a  major  factor  that  complicates  the  obtaining  of  a  solu- 
n,  l 

tion  to  the  problem. 

The  term  p_  .  is  the  probability  that  no  machines  are 
U ,  1 

down  in  time  period  i  .  Since  this  term  is  a  common  factor  in 

all  expressions  of  p  .  ,  it  acts  as  a  normalizing  factor  and 

n,i 

assures  that  the  sura  of  all  probabilities  over  all  states  equals 
one. 


Equation  (II-2)  gives  the  mean  failure  rate  A^  for  two 

cases,  one  when  the  number  of  machines  in  the  system  is  increasing, 
and  one  when  they  are  decreasing.  When  the  population  is  increas¬ 
ing  (M^  >  )  ,  the  expression  is  composed  of  three  terms.  The 

first  reflects  the  failure  rate  of  all  components  introduced  in 
time  period  i  ,  and  the  second  term  reflects  the  change  in  failure 
rate  due  to  those  components  that  failed  and  were  repaired  in  the 
previous  time  period  (i-1)  .  The  last  term  is  associated  with 
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chose  components  that  were  neither  introduced  in  time  period  i  nor 
repaired  in  time  period  i-1  ;  such  components  maintain  the  previous 
period's  failure  rate.  It  is  assumed  that  components  that  are 
either  introduced  or  repaired  in  a  time  period  i  attain  the  same 
failure  rate.  The  expression  for  the  mean  failure  rate  when  the 
population  is  decreasing  (M^  <  is  similar  to  the  expression 

for  an  increasing  population  except  for  the  first  term.  Since  there 
are  no  machines  being  introduced  in  time  period  i  ,  there  is  no 
dependence  on  the  reliability  of  new  machines. 

The  use  of  a  weighted  mean  for  the  expected  failure  rate  was 
investigated  in  References  [4]  and  [5].  It  was  found  that  this 
yields  conservative  results  when  the  availability  is  calculated. 

The  validity  of  the  approximation  decreases  as  the  failure  rate  in¬ 
creases.  In  the  use  of  the  approximation,  it  is  assumed  that  no 
machine  fails  more  often  than  once  in  a  time  period. 

Equation  (II-3)  gives  the  mean  number  of  machines  repaired  • 
in  a  time  period. 

Next,  the  constraint  on  the  availability  of  spares  must  be 
specified.  This  requires  the  calculation  of  the  failure  point 
probabilities,  q^  ,  which  are  the  state  probabilities  determined 

at  the  instant  of  failure.  The  state  probability,  p^  ,  is  the 

probability  that  n  machines  are  down  at  any  arbitrary  time.  The 
q  ' s  ,  on  the  other  hand,  are  measured  only  at  the  instant  of  fail¬ 
ure.  Between  failures  it  is  possible  for  a  machine  to  be  repaired, 
thereby  changing  the  state  of  the  system.  The  failure  point  prob¬ 
abilities,  therefore,  are  the  state  probabilities,  conditioned  on 
the  instant  of  failure;  or,  using  the  standard  terminology  of  the 
literature, 

qn  =  Pr[n  in  system |a  failure  is  about  to  occur]  . 


v  ■  ....  •  v  r.  *  •■■•v-  v  ..  ...  *  '  ft  .A  »» ft.  atAiyat,  -V* 
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Since  Che  constraint  on  the  availability  of  spares  requires  that 
a  spare  gas  turbine  be  available  when  requested  with  a  probability 
of  0.90,  and  that  the  request  be  made  at  the  time  of  failure,  the 
availability  constraints  are  formulated  in  terms  of  the  q^'s  • 

From  Bayes'  theorem, 

Prfn  in  system]  Pr[ failure  about  to  occur jn  in  system! 

^n  £  [Pr{n  in  system}  Pr{failure  about  to  occur|n  in  system}] 
n 


For  this  problem,  q^  ^  is  determined  to  be 


i  n,  i 


*i+Mi 


(0  <n  <  y .) 

— i  l 


mi  ■ .?  u'yi>p*.i 


*n,i 


(M1-iHy1)pn<1 


(II-4) 


yi+Mi 


,  (y<n<yi+Mi) 


m  -  £  i 

x  o=«  *  X'*1 


£=y 


where  the  p^  ^  are  given  in  Equation  (II-l).  Therefore,  the 
availability  constraint  is  specified  as 


y.  q  ,  >  0.90  ,  i=l,2, 

JTn  n> 1  — 


>k 


where 


q  ^  =  the  steady  state  probability 
n*  that  n  machines  are  down  at 
the  time  of  a  machine  failure, 
year  i  ;  and 

k  =  length  of  the  planning  horizon. 


II. 3  Cost  Model 


The  problem  as  originally  stated  in  Reference  [1]  had  a  cost 
model  that  was  dependent  on  four  types  of  costs:  purchase  cost  of 
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spares,  purchase  cost  of  repair  channels,  repair  costs,  and  invest¬ 
ment  costs  for  component  improvement.  This  cost  function  is  ex¬ 
pressed  as: 


where 

(a)+  =  max(a,0) 

Cf  i  =  purchase  cost  of  a  repair  channel,  year  i 

C9  .  =  purchase  cost  of  a  spare  component,  year  i 
z ,  1 

C.  .  =  cost  of  repairing  a  component,  year  i 

j,  i 

^  =  investment  in  reliability  growth  (CIP)  pro- 
'  gram,  year  i 

d  =  discount  factor 

k  »  number  of  years  in  planning  horizon. 

The  mathematical  programming  problem  becomes 


c j , c2 , . 


cM1"  y  ,  z '  &  di[ci,i<crVi>++c2,i(yr>'i-/ 

’ V  yl*y2'm  *yk  11 

+  C  .R.  +C.  . 

3,i  l  4,iJ 


V1 


Subject  to  53  q  >0.90  ,  (i=l,2, . . . ,k) 

n=0  0,1 


c^  >  0  ;  integer 
>  0  ;  integer. 


C II— 5) 


(II-6) 


Since  the  last  term. 


is  constant  for  each  time  peri¬ 


od,  it  can  be  dropped  from  the  optimization.  The  third  term,  which 
is  dependent  on  ,  is  an  explicit  function  of  y^  and  an  im¬ 
plicit  function  of  c^  and  y^  .  In  the  original  analysis  of 
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Reference  [1]  this  term  is  dropped  from  the  cost  function.  Since 
R  is  relatively  constant  [if  (II-6)  is  satisfied],  this  term  will 
~nly  have  a  secondary  effect  on  optimality.  The  optimization 
problem  then  becomes 


Min 

vyi 


+  C2,i(!'i-yi-l)+] 


V1 

subject  to  V  q  .  >  0.90  ,  for  all  i 
n=0  n>1 

c  >  0  ;  integer 

y  >  0  ;  integer. 


Since  there  is  no  cost  benefit  associated  with  selling  a 
spare  or  repair  channel,  the  problem  can  be  reformulated  to  have  a 
linear  objective  function,  which  may  be  written  as  follows: 


cl,c2,-. 


Min 

V  yi’y2 


k 

E 

i=l 


subject  to  Y'  q  .  >  0.90  ,  for  all  i 
n=0  0,1 


c^  integer 

c  .  >  c .  , 
x  =  l-l 

Ci  >  0 

yi  integer 
yi  =  yi-l 


y1  >  o  . 


(n-7) 


This  is  the  basic  programming  problem  whose  solution  is  dis¬ 
cussed  in  this  paper.  It  is  a  nonlinear  integer  programming  problem 
with  a  linear  objective  function.  A  solution  is  achieved  here  using 


15 


an  amalgamation  of  partial  enumeration  and  branch  and  bound,  and 
dynamic  programming  techniques.  The  techniques  developed  can  be 
applied  to  a  more  general  objective  function  as  well  as  to  a  more 
general  problem  type,  which  will  be  discussed  in  the  final  chapter. 

II. 4  Heuristic  Approach 

An  attempt  to  solve  this  problem  using  a  heuristic  approach 
was  presented  in  [1].  This  approach  is  predicated  on  a  procedure 
that  minimizes  the  incremental  cost  of  purchasing  spares  or  repair 
channels  for  each  successive  time  period.  Starting  with  the  number 
of  spares  and  repair  channels  already  purchased  in  the  previous  time 
period,  the  algorithm  searches  for  the  policy  that  satisfies  the 
availability  constraint  for  the  current  time  period  for  a  minimal 
additional  cost.  If  in  a  time  period  the  availability  constraint 
is  satisfied  by  the  policy  of  the  previous  time  period,  a  reduction 
procedure  is  employed  that  tries  to  decrease  the  number  of  spares 
or  repair  channels  in  such  a  manner  that  the  availability  constraint 
is  merely  satisfied,  even  though  there  is  no  explicit  cost  savings. 
External  manipulation  of  the  input  to  the  algorithm  allows  modifica¬ 
tion  of  the  policy  used  in  previous  time  periods  based  upon  the  poli¬ 
cy  obtained  for  later  time  periods. 

Based  upon  a  comparison  of  computational  experience,  this  al¬ 
gorithm  provides  results  that  approach  the  optimal  policy,  but  does 
not  assure  that  the  optimal  policy  will  be  obtained.  Its  basic 
weakness  is  its  one  year  at  a  time  "optimality,"  which  does  not 
foresee  future  cost  savings.  The  computational  results  are  pre¬ 
sented  in  Chapter  IV  of  this  dissertation. 


CHAPTER  III 


METHODOLOGY 

III.l  Introduction 

In  this  chapter  the  methodology  developed  to  solve  this 
spares  provisioning  problem  is  presented.  There  are  many  unique 
aspects  of  the  problem  that  require  individual  solution  techniques; 
the  difficulties  created  by  these  aspects  and  the  techniques  re¬ 
quired  to  overcome  them  are  discussed  here.  The  overall  solution 
technique  and  proof  of  the  validity  of  the  method  are  presented 
as  well. 

Contained  in  this  chapter  are:  a  discussion  of  each  of  the 
problematic  areas,  an  overview  of  the  solution  technique,  the  method 
for  generating  regions  that  satisfy  the  availability  constraint  for 
each  time  period,  the  method  of  combining  these  regions  across  time 
periods  to  form  the  feasible  space,  the  dynamic  programming  solution, 
and  a  discussion  of  the  branch  and  bound  algorithm  utilized. 

III. 2  Problematic  Areas 

The  most  obvious  difficulty  of  this  problem  is  that  it  is 
a  nonlinear  integer  program.  A  second  difficulty  is  associated  with 
the  complicated  functional  dependence  in  the  equations  defining  the 
problem.  The  last  major  complication  is  that  a  key  parameter  cannot 
be  predicted  until  the  optimal  policy  is  specified. 

The  nonlinear  integer  aspect  of  the  problem  is  obvious. 

While  the  objective  function  can  be  restructured  into  a  linear  form 
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by  adding  Che  proper  constraints,  the  nonlinearity  of  the  state 

p  affects  the  determination  of  the  constraint  equations.  As 
n.i 

will  be  discussed  later,  the  linearity  of  the  objective  function 
is  important  in  the  solution  technique. 

The  form  of  the  equation  of  state  as  seen  in  Equation  (II-l) 
requires  an  enumerative  approach  to  the  solution.  Unless  the  values 
of  c  ,  the  number  of  repair  channels,  and  y  ,  the  number  of  spares, 
are  known  relative  to  M  (the  number  of  turbines  on  line)  and  n 
(the  number  of  turbines  down),  it  is  not  possible  to  specify  the 
form  or  value  of  p^  ,  the  probability  that  n  turbines  are  down. 

However,  for  any  given  "policy"  (c^,^, — ,  y^.y,,,...)  ,  it  is  pos¬ 
sible  to  evaluate  p  .  for  all  n  . 

1  «  -» 


The  remaining  difficulty  arises  from  the  nature  of  A ,  the 

average  failure  rate  in  time  period  i  .  As  used  in  Equation  (II-l), 
it  appears  that  Ai  is  a  parameter  that  is  required  to  determine 


the  state  variables  p 


.  .  Actually,  A.  is  an  implicit  function 
n,  1  i 


of 


and 


yj 


for  j  less  than  i 


With  A^  dependent,  bound¬ 


ing  the  feasible  region  in  time  period  i  becomes  extremely  diffi¬ 
cult,  if  not  impossible.  In  using  a  branch  and  bound  procedure, 

A^  is  approximated  in  a  manner  that  makes  it  independent  of  the 


^Ci’^0  choices  of  the  previous  time  periods.  By  employing  such 

a  procedure  it  is  possible  to  specify  the  convex  hull  of  the  region 
that  satisfies  the  availability  constraints  for  each  time  period 
independent  of  the  other  time  periods. 


Critical  to  the  approach  used  for  solving  this  problem  is 
an  understanding  of  the  nature  of  the  solution  to  problems  with 
linear  objective  functions.  Theorem  2  of  Reference  [6]  states  that: 

"...  a  linear  objective  function  assumes  its 
minimum  at  an  extreme  point  of  the  convex  set 


K  generated  by  the  set  of  possible  solutions 
to  che  linear  programming  problem.  If  it  as¬ 
sumes  its  minimum  at  more  than  one  extreme 
point,  then  it  takes  on  the  same  value  for 
every  convex  combination  of  these  particular 
points . " 

An  extreme  point  of  a  convex  set  is  a  point  in  that  set  that  cannot 
be  expressed  as  a  proper  convex  combination  of  any  two  distinct 
points  in  the  set.  The  fact  that  this  problem  is  not  a  linear  pro¬ 
gramming  problem  in  the  normal  sense  does  not  invalidate  the  theorem 
above.  The  convex  hull  of  the  lattice  of  integer  pairs  satisfying 
the  availability  constraint  could  be  determined  and  specified  by  a 
finite  number  of  linear  inequalities.  Since  the  remaining  con¬ 
straints  (ci=ci_i’  yi^yi.-i’  c  >0,  are  linear>  the 

problem  could  then  be  formulated  as  a  linear  integer  programming 
problem. 

By  the  above  theorem,  however,  it  is  not  necessary  to  go 
through  the  cumbersome  technique  of  solving  the  linear  programming 
problem.  In  the  approach  described  herein,  the  extreme  points  of 
the  convex  hull  of  the  lattice  of  integer  pairs  satisfying  the 
availability  constraint  are  determined  explicitly  through  the 
application  of  a  specified  algorithm.  Once  these  extreme  points 
have  been  enumerated,  the  complete  set  of  extreme  points  associated 
with  the  remaining  linear  constraints  linking  the  time  periods  can 
also  be  listed.  A  dynamic  programming  approach  is  then  used  to 
determine  the  extreme  point  that  gives  the  minimum  of  the  objec¬ 
tive  function. 

Therefore,  the  overall  approach  can  be  considered  in  terms 
of  the  three  major  techniques  of  which  it  is  composed: 

(1)  a  branch  and  bound  procedure  based  upon  the  approx¬ 
imation  of  A^  .  By  approximating  A^  ,  the  de¬ 
pendence  of  A^  on  the  c  j ' s  anc*  ^j*S  ’  ^or  ^ 
loss  than  i  ,  can  be  eliminated  from  the 
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subproblems  of  the  branch  and  bound  technique. 

•  _ 

In  this  manner,  X^  becomes  a  parameter  of 
the  problem. 

(2)  A  finite  enumeration  technique,  made  possible 
by  a  linear  objective  function,  which  allows 
the  solution  to  lie  on  the  extreme  points  of 
the  convex  hull.  These  extreme  points  are  de¬ 
termined  explicitly. 

(3)  Dynamic  programming  techniques,  applied  to  ob¬ 
tain  an  optimum  solution  for  the  given  set  of 
X/s  ,  once  the  set  of  extreme  points  of  the 

convex  hull  has  been  determined. 

The  remainder  of  this  chapter  will  be  devoted  to  the  enumeration 
technique  and  the  proof  that  the  procedure  is  finite,  review  the 
dynamic  programming  approach,  and  show  that  branch  and  bound  is 
applicable  to  this  problem.  The  branch  and  bound  procedure  is 
discussed  next. 


III. 3  Branch  and  Bound 


In  determining  if  the  availability  constraint  is  satisfied 
at  a  particular  pair,  it  is  necessary  to  evaluate  X^  . 

Every  p  .  ,  the  state  probability  for  time  period  i  ,  is  depend- 
n,  l 

ent  on  X^  .  It  would  be  convenient  if  the  X^'s  were  parameters 

of  the  problem,  but  wc  know  that  they  are  dependent  on  the  choice  of 
policy  in  the  previous  time  periods.  The  mean  failure  rate  for  a 
time  period  i  is  given  by  Equation  (II-2)  and  repeated  here: 


\  ■ 


l'Wi>Ai  +  “t-ih-i  +  •  Wi 

I'i-i'm  +  “i-rVi’Vil/Vi  ■ 


As  was  mentioned  in  Chapter  II,  this  is  an  approximation  to  the 
failure  rate  associated  with  the  process.  Since  the  approximation 
deteriorates  as  A^  increases,  it  has  therefore  been  assumed  that 

is  small  and  that  thus  A^  is  also.  An  upper  limit  imposed  on 

A  is  given  later  in  this  section. 

The  dependence  of  A^  on  the  policies  of  previous  time 

periods  is  introduced  through  the  ,  the  expected  number  of • 

machines  repaired  in  the  previous  time  period.  Through  its  de¬ 
pendence  on  p  .  ,  R.  is  both  explicitly  dependent  on  y.  and 

implicitly  dependent  on  y_^  and  c^  .  If  R^  were  a  parameter 

of  the  problem,  A^  ,  for  all  i  ,  could  be  specified.  The  branch 

and  bound  procedure  is  based  upon  the  parameterization  of  R^  . 

The  branch  and  bound  procedure  is  predicated  on  the  assump¬ 
tion  that  the  initial  problem  may  be  decomposed  into  a  set  of  sub¬ 
problems  with  a  subset  of  these  problems  solved  in  a  series  of 
s  tages  [ 7 ] . 

Suppose  the  original  problem  is  of  the  form 

rain  f(x) 
s.  t.  x  £  S  , 

and  X*  is  the  solution  set.  A  branch  and  bound  procedure  ini¬ 
tially  addresses  an  extended  problem  of  the  form 

min  F(x) 

s. t.  x  £  T  , 

where  S  C  T  and  f(x)  =  F(x)  for  x£  S  .  A  sequence  of  stages 
is  created,  the  first  stage  defined  by  the  initial  extended  problem 
above.  Associated  with  a  subsequent  stage  n  is  a  collection  of 
problems  of  the  form 
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min  F(x) 

s.t.  xeT1  ,  i  e  N(n)  . 

Here  N(n)  defines  the  set  of  indices  for  the  collection  of  prob¬ 
lems  composing  stage  n  . 

The  branch  and  bound  method  will  ensure  that 

U  T1  C  T  and  X*  H  (  U  T*\  t  $ 
i£N(n)  \icN(n)  / 

at  each  stage  n  ;  i.e.,  each  stage  is  comprised  of  subproblems 
whose  collective  feasible  region  contains  at  least  one  solution  of 
the  original  problem.  The  solution  to  each  of  the  i  problems 

provides  a  lower  bound  (LBK  to  f  over  SHT1  .  The  best  lower 

bound  at  stage  n  ,  (BLB)  ,  is  the  minimum  of  the  set  of  lower 

n 

bounds  at  stage  n  ,  written  as 

(BLB)  =  min  (LB).  . 

n  ieN(n)  1 

Since  SCT  and  f(x)  *  F(x)  for  xES  ,  the  (BLB)  <  f(x*)  , 

~  n  - 

with  x* E  X*  .  At  each  stage  of  the  branch  and  bound  procedure  a 

new  set  of  problems  is  solved.  This  set  of  problems  is  chosen  in 

a  manner  such  that  * 

U  T1  C  U  T1  . 

i£N(n+l)  ieN(n) 

Proceeding  in  this  manner  from  stage  to  stage,  successive  best  lower 
bounds  arc  generated,  leading  to  the  following  property  of  the  branch 
and  bound  procedure: 

(BLB),  <  ...  <  (BLB)  <  f(x*)  . 

i  ~  n 


IF  at  any  stage  of  the  problem  a  solution  that  gives  the  best 
lower  bound  at  that  stage  is  contained  in  S  ,  the  feasible  region 
of  the  original  problem,  then  this  solution  is  also  a  solution 
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to  the  original  problem.  Also,  a  lowest  valued  solution  that  is 
feasible  to  the  original  problem  gives  the  best  upper  bound  (BUB) . 

If  the  BUB  =  BLB  ,  the  problem  is  solved  with  the  solution  value 
equal  to  the  BUB  =  BLB  . 

In  applying  the  branch  and  bound  procedure  to  the  spares 
provisioning  problem,  we  use  the  objective  function  as  given  in 
(II-7).  The  feasible  region  of  the  extended  problem  is  created  by 

-  A 

modifying  the  definition  of  X^  .  By  defining  X^  such  that 

X^  <  X^  for  all  feasible  X^  and  for  all  i  ,  we  generate  a  region 

T  such  that  SCT  . 


/v  - — 

Theorem  1:  For  X^  <_  X^  ,  a  feasible  region  T  associated  with 

the  problem  given  in  (II-6)  can  be  generated  such  that  SCT  , 
where  S  is  the  feasible  region  of  problem  (II-6). 


Proof:  This  can  be  proven  by  considering  the  availability 

constraint.  This  is  the  only  constraint  that  is  affected  by  modi¬ 
fying  X^  .  Rewriting  the  availability  constraint  as  in  Appendix 

B,  Equation  (B-l) ,  we  have 


From  (II-l),  p^  ^  is  proportional  to  A? 
fore,  the  following  notation  is  introduced: 


for  all  cases. 


There- 


P 


n,  i 


=  Q  .XV 

n,  1  l 


* 


where  Q  is  a  function  of  M,  ,  u.  ,  c,  ,  and  y.  .  In  that 
n,i  i  i  i  i 

the  time  period  remains  the  same  for  the  duration  of  this  discus¬ 
sion,  the  subscript  i  will  be  dropped  to  ease  the  reading.  The 
availability  constraint  can  be  rewritten  as 


A 
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M-l 

E 

l-a  1-0  _ 

a  -  y-1 

o  \° 

n»0 

Assume  that  the  availability  constraint  is  satisfied  for  some 
c',y*  .  We  will  now  show  that  the  pair  (c’,y')  also  satisfies 


l-a 


lo  “‘H,  ^ 


>  21 

a  —  y-1 

E  Qnxn 

n=0  n 


provided  that  A  <  A  .  Equivalently,  it  will  be  shown  that 

M  i  .  .  M  ■ 

r  A^  £  (!-£)<>  Xy+J 

i+y  £n  Mj+y 


j=0 


y-i  __ 

£  Qn^ 


>  A* 


n-0 


£  0-r,  X 


n=0 


or 


(I  w-*v*i  jg  vi  i  jg  <^<wi  jg 


M  y-1 


M  y-1 


£  £  £  £  <l-s*WVy+J*"  i  • 


j=0  n=0 


j=0  n=0 


4  -¥)<■• 


but  for  A/A  >  1  and  y+j  >  n  , 

n-\y+J 


r-dj 


>  0  . 


Since  all  other  terms  in  the  summation  are  greater  than  0,  this 
inequality  holds.  Therefore,  if  the  inequality  associated  with 
the  availability  constraint  holds  for  a  pair  (c',y')  and  \  , 
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then  it  also  holds  for  the  same  pair  (c',y')  when  X  <  X  .  This 
implies  that  for  any  time  period,  any  point  that  satisfies  the 

availability  constraint  using  the  actual  X^  will  also  satisfy  the 

^  _  « 

availability  constraint  using  X^  <  X^  so  if  the  extended  problem 

is  associated  with  using  a  suitable  set  of  X^  we  have  SCT  . 

Here  S  is  the  feasible  region  of  the  original  problem,  using  X^ 

for  evaluating  the  availability  constraints,  and  T  is  the  feasible 
region  associated  with  the  extended  problem,  using  X  <  X^  for 

evaluating  the  availability  constraints. 


The  problem  remains  to  determine  a  suitable  set  of  X^ 
values  .  From  Equation  (II-2)  we  see  that  X^  and  X 2  are  inde¬ 
pendent  of  any  policy  and  are  dependent  only  on  the  specified  pa¬ 
rameters.  Thus, 


X1-X1 


ix2  +  (M1/M2)(X1-  X2)  ,  M2  >  Hj_ 

X1  *  M2  <  Mi  ' 

For  subsequent  time  periods  X^  is  specified  by: 

Xl  \\-i +  (Vi'ViHVi-Vii  •  Wi  ’ 


R  -  X 


r  Vyi  i 

b '  4i«  (nr,i>,-iJ 


Introducing  the  notation 


mj^i 

<1  -  E  ("-Ri'Pn.i  ■ 

n^y^l  ’ 
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we  get 


(Ai+<Mi-i/"i)ai-rV +  xi-iai-rxi-i)<Mi-rxi-i)  '  \  •  Wi 
1 1  (xi-i  +  xi-i(Xi-rxi-i)(Mi-rxi-i>  1  Mi-i  ■  Wi 


(III-l) 


We  now  have  A^  totally  dependent  on  the  choice  of  the  policy  of 
the  previous  time  periods  contained  in  the  ^  .  Note  that  A^ 
is  not  dependent  on  .  The  bounds  on  A^  can  be  determined. 


Limits  on  X.  can  be  determined  by  noting  that,  since  p 

l  n ,  l 

M,+y. 

is  a  probability,  it  follows  that  0  <  p  .  < 1  and  )  n  p  .  *  1  • 

n ,  i  —  n=u  n ,  i 

Also,  in  the  summation  defining  X^  ,  we  have  y ^+1  n <  IL+y  . 

M.+y, 

Therefore,  for  X.  =  (n-y^.  , 

Mi+yi 

0  <  X  <  M,  £  p  . 

-  i  —  i  ,,  n,i 

n=yi+l 

or 

0  Jx.  <M.  . 

To  determine  the  initial  extended  problem  in  the  branch  and 
bound  procedure,  we  must  find  a  A^<  A^  »  therefore,  we  will  choose 

A  _ 

a  value  of  X.  ,  such  that  A.  =  min  (A.)  .  Since  0  <  X .  ,<11.  ,  , 
l-l  i  i  *=  i-1  —  i-I 

Xi-1 

from  (III-l)  the  rule  for  choosing  A^  should  be  determined  either 
by  setting 


xi-i  '  0  • 


X.  ,  <  A.  , 

1-1  ^  1-1 


X.  .  =  M,  ,  A.  >  X.  . 
i-I  i-1  i-I  i-I 
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However,  the  problem  is  more  complex  than  that.  All  sub¬ 
sequent  Aj's  ,  j  greater  than  i  ,  are  dependent  on  the  value  of 

.  While  may  be  chosen  to  minimize  A^+^  *  it  “ay  maximize 

\^+2  >  an<*  so  on* 

By  using  the  explicit  dependence  of  A^  on  A^  we  can 

specify  a  condition  on  A^  to  ensure  that  the  minimization  of  A^ 

will  also  minimize  A^  ,  for  j  greater  than  i  .  Since  A^  is 

a  parameter  of  the  problem,  if,  when  the  problem  is  specified,  the 
A^'s  do  not  satisfy  this  bound,  there  is  no  assurance  that  the 

procedures  described  herein  will  work. 


First  note,  from  expression  (II-2) ,  that  A^  is  an  average, 

and  that,  therefore,  A  £  max  (A.)  .  Next  we  examine  the  behavior  of 

1  2 

A^  •  If  A^  is  a  monotonically  increasing  function  of  A^_^  ,  then 

we  know  that  the  minimum  value  of  A.  will  occur  at  the  minimum 

l 

value  of  1.  ^  .  Examining  (III-l) ,  we  see  that  Ai  is  quadratic 

in  A.  .  if  X.  <  M  ,  and  linear  in  A.  ,  if  X.  ,  =  M.  .  . 
i-l  l-l  i-l  l-l  l-l  l-l 

Therefore,  for  ^  ,  since  we  have  monotonicity ,  the  minimum 


A.  is  attained  when 
i 

A i-l  ls 

minimum. 

For 

X.  <  M.  .  ,  if  A. 
i-l  i-l  i-l 

is  constrained  to  be 

less  than 

*i*-l  • 

where 

3Ai(A*_i)/3A*-i  *  0  , 

then  the  minimum  X^ 

is  attained  when 

Xi-1 

is  minimum;  i.e., 

l(li  f)  is  a  monotonically  increasing  function  of  A^_^  for 


Evaluating  3A^/3A^_^  using  (III-l), 


27 


ax. 


M.  i  M  -  X 


Mi-1 


M  -  X 

1  +  \_L~  wi-i '  2Xi-i)  • 


Mt  <  Mi-i 


Setting  aX./aX.  to  zero  to  determine  X*_^  »  we  Set 

\l  ’  *  V  •  • 


Since 


0  =  Xi-1  <  Mi-1  ’ 


,  X.  . 

xi-i  ^i  +  -yi  • 


then 


Therefore,  if 


Xt  <  y  ,  for  all  i  ,  since  X^  <  X^  , 

then  we  know  that  X^<y  for  all  i  .  Therefore  X^_^  can  never 
be  greater  than  X* ,  so  X^  is  a  monotonically  increasing  func¬ 
tion  of  X.  ,  and  the  minimum  X.  will  be  attained  for  the  mini- 
l-l  l 


mum  X .  , 
l-l 


This  is  a  sufficiency  condition.  Tighter  requirements  could 
be  found  for  the  upper  bound  of  X_^  ^  to  assure  that  the  minimum 

value  of  X^  is  attained  at  the  minimum  value  of  ^  j  .  In  that 

the  expression  (II-2)  is  an  approximation  to  the  actual  failure  rate 
of  the  process,  and  the  approximation  is  good  only  for  small  values 
of  X^  ,  this  sufficiency  condition  does  not  impose  any  additional 

practical  restrictions  on  the  problem.  In  the  next  chapter  we  will 
see  that  the  numerical  data  more  than  satisfy  this  condition. 
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It  is  now  possible  to  specify  the  branch  and  bound  proce¬ 
dure.  The  original  problem  is  • 


k  .  _  _ 

Min  £  d1  C.  .c.  +  C  .y 

:.,y.,vi  i=l  L  1*1  1  2’x  iJ 

1  1 


yi_1 


t.  53  q  .  ^  0.90  ,  'v'i 


n=0  0,1 


c  .  >  c .  , 

l  -  i-l 


yi  =  yi-l  * 


'■57'i 

N/i 


y1  >  o 


ci  >  o 


c^  and  y^  are  integers  for  all  i  . 


The  q  .  are  implicit  functions  of  X.  which  are  them- 
n,  l  l 

selves  implicit  functions  of  c^ ,y^  for  j  less  than  i  .  The 

-  A  a 

extended  problem  is  created  by  replacing  X^  by  X^  .  The  X^  are 

A  — —  A 

independent  of  c^.y^  and  X^<X^  "  '^1e  \  are  8iven  by: 

( III-2a) 

xi +  <Mi-i/Mi) <xi-rxi> +  Si-i<Mi-i/V <xi-r Vi>  > 


for  X.  ,  >  X.  ,  , 
i-l  i-l 


Ai-1+Ai-i(Ai-rxi-i)  * 


Mi  <  Vi 


or  for  X.  ,  <  X.  ,  , 
i-l  i-l 


X.  = 
i 


\l  +  (Mi-i/Mi)(Ai-rV  •  Mi^Mi-i 


( III-2b) 


'i-l  ’ 


M  <  M  , 
i  i-l 


I 
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It  has  been  shown  in  this  section  that  a  valid  extended  problem  is 
crea-  2d  by  this  substitution. 

The  extended  problem  is  solved  by  techniques  which  will  be 
explained  in  the  following  sections.  If  the  solution  is  feasible, 
the  problem  is  solved.  If  the  solution  is  not  feasible,  branching 
occurs . 


If  the  set  {(Cj,y^);  j  < i}  is  known,  it  is  possible  to 

calculate  exactly.  Therefore,  if  a  set  {(c^,y^);  j  < i}  is 

specified  for  all  time  periods,  the  A^  for  all  time  periods  can 

be  calculated  exactly,  and  this  policy  can  be  tested  to  determine 
if  it  is  feasible. 


The  branch  and  bound  technique  when  applied  to  this  problem 
proceeds  in  the  following  manner.  The  first  subproblem  is  generated 

A 

by  calculating  A^  for  all  time  periods.  For  time  periods  1  and  2, 
A^  *  Af  ,  since  for  these  time  periods,  A^  is  dependent  only  on 
the  input  parameters  and  is,  therefore,  known  exactly.  For  all 

A 

other  time  periods  A^  is  calculated  using  Equations  (III-2a)  or 

(III-2b).  This  subproblem  is  solved  using  the  techniques  described 
in  the  following  sections.  If  the  solution  to  this  subproblem  is 
feasible,  we  have  found  the  solution  to  the  original  problem  and  we 
are  finished.  If  the  solution  is  infeasible,  branching  is  performed. 

In  solving  the  first  subproblem,  a  set  of  m  pairs 
{(^2*^2^?.’  is  determined  for  the  second  time  period.  We 

know  that  the  solution  to  the  initial  problem  must  contain  one  of 
these  pairs  as  the  optimal  pair  for  the  second  time  period.  This 
is  because,  in  the  procedure  that  is  described  in  the  following 
sections,  for  a  given  A^  all  of  the  extreme  points  of  the  region 

satisfying  the  constraints  are  enumerated.  This  enumeration  is 

/A 

dependent  on  A^  ,  but  for  the  second  time  period  (as  well  as  for 

/\  _ 

the  first  time  period)  ^2  =  ^2  ' 


If  there  are  m  such  pairs,  then  for  each  one  of  these 

—SL 

pairs  {(c^,y^) V  »  a  A3  can  be  calculated.  If  the  optimal 

policy  associated  with  the  solution  of  the  initial  problem  contains 

_2  _ 

<c2’^2^ l  *  then  ^3  is  the  exact  ^3  for  the  third  time  period. 

— £ 

Thus  we  have  created  m  new  subproblems,  each  having  a  A^  asso¬ 
ciated  with  each  of  the  m  pairs  ^c2,y2^Z  '  For  a11  subsequent 
time  periods  (>4)  for  each  subproblem,  the  ^  are  calculated 
using  (III-2a)  or  (III- 2b). 

Each  new  subproblem  is  solved  using  the  solution  techniques. 
The  problem  that  has  the  best  lower  bound  is  then  tested  for  feasi¬ 
bility.  If  it  is  feasible,  the  solution  has  been  found.  If  it  is 
not  feasible,  branching  continues  based  upon  the  enumerated  set  of 
pairs  ^3*^3)  2  '  ^  the  solution  to  the  initial  problem  contains 

the  (c,y)  pair  of  the  second  time  period  that  generated  the  sub- 
problem  from  which  the  current  branching  is  to  occur,  then  the  op¬ 
timal  policy  for  the  third  time  period  for  the  solution  of  the  ini¬ 
tial  problem  must  come. from  the  set  of  pairs  (c„,y  ).  . 

J  J 

The  new  branches  are  generated  in  the  same  manner  previously 
—2 

performed,  with  an  exact  A^  calculated  for  each  of  the  pairs 

^C3’^3^£  ’  Tbcre  is  a  new  subproblem  associated  with  each  of  the 
— 2 

A ^  ,  with  the  A^  calculated  for  subsequent  time  periods  (^5) 

using  \III-2a)  and  (III-2b).  After  the  branching  a  new  best  lower 
bound  is  determined  and  tested  for  feasibility. 

This  procedure  continues  until  a  feasible  best  lower  bound 
is  found.  Branching  is  always  based  upon  the  set  of  pairs  (e^.y 

for  the  last  time  period  in  the  subproblera  associated  with  the  best 

lower  bound,  for  which  the  exact  A.  has  been  used.  The 

1 
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_  i 

determination  of  an  exact  is  dependent  on  the  specification  of 

(c,y)  pairs  for  all  previous  time  periods. 

This  procedure  is  equivalent  to  taking  all  possible  combina- 
tions  of  policies  for  each  time  period  based  upon  the  enumeration 
^of  those  points  that  satisfy  the  constraints.  Since  it  is  not  pos¬ 
sible  to  enumerate  correctly  the  set  of  policies  in  a  given  time 
period  until  the  policies  of  all  the  previous  time  periods  have 
been  specified  (A^  cannot  be  calculated),  it  is  necessary  to  do 

this  enumeration  in  a  sequential  order.  The  advantage  of  the 
branch  and  bound  procedure  is  that  it  offers  the  hope  that  the  so¬ 
lution  to  the  initial  problem  will  be  found  before  it  is  necessary 
to  go  through  all  of  the  branches.  This  hope  has  been  justified 
in  practice. 

Since,  for  a  given  A^  ,  we  are  interested  in  only  a  finite 
number  of  pairs  (c^.y^)  that  satisfy  the  constraints  of  the  prob¬ 
lem,  as  will  be  discussed  in  the  next  section,  and  since  there  is 
a  finite  number  of  time  periods,  there  can  only  be  a  finite  number 
of  branches.  Therefore,  the  branch  and  bound  procedure  must 
converge  [8]. 

The  procedure  can  now  be  summarized  as  follows: 

(1)  Create  the  extended  problem  by  replacing  X^ 

A 

by  the  approximation  A^  ,  as  given  in  Equa¬ 
tions  (III-2a)  and  (III-2b). 

(2)  Solve  the  extended  problem. 

(3)  If  the  solution  is  feasible,  the  problem  is 
solved. 

(4)  If  the  solution  is  infeasible,  branch. 

(5)  Based  upon  the  last  time  period  &  for 

A  - 

which  A  »  A.  for  all  k  less  than  £  , 
i  k  * 
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calculate  the  set  .  Create  a  new  set  of 
‘  problems. 

(6)  Solve  the  new  set  of  problems. 

(7)  If  the  problem  associated  with  the  BLB  is 
feasible,  it  is  the  exact  solution  and  the 
problem  is  complete. 

(8)  If  the  problem  is  infeasible,  branch  as  in 
Step  5. 

We  will  next  look  at  the  techniques  that  are  required  to 
solve  the  extended  problem. 

III. 4  Enumeration  Techniques 

The  objective  function  has  been  manipulated  into  a  linear 
form;  therefore,  if  there  is  a  bounded  solution  then  there  is  one 
which  will  be  on  the  convex  hull  of  the  feasible  region.  For  the 
branch  and  bound  approach  that  is  being  used,  it  Is  also  true  that 
for  any  stage  of  the  procedure  the  solution  will  be  found  on  the 
convex  hull  of  the  feasible  region  of  a  subproblem.  A  method  of 
generating  the  convex  hull  is  therefore  required.  The  enumeration 
technique  is  the  first  step  of  that  process. 

For  each  time  period,  and  for  each  subproblem  of  the  branch 
and  bound  method,  the  convex  hull  of  the  region  satisfying  the 
availability  constraint  can  be  enumerated.  In  this  section  it  will 
be  shown  that  the  enumeration  is  finite.  Also  presented  in  this 
section  is  the  algorithm  used  in  the  enumeration.  Because  the  fol¬ 
lowing  discussion  is  independent  of  time,  the  subscript  i  will 
again  be  dropped  to  simplify  notation. 

Figure  III-l  illustrates  the  procedure  used  to  enumerate 
the  region  satisfying  the  availability  constraint.  There  are  three 
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lines  of  interest  on  this  figure:  the  solid  step-like  line,  the 
dashed  line,  and  the  diagonal  line  c=y  .  Each  of  these  is  sig¬ 
nificant  in  the  discussion  of  the  enumeration  process. 

Because  of  the  physical  reality  of  the  problem  (fractional 
spares  or  repair  channels  are  meaningless)  and  the  nature  of  the 
equation  defining  the  state  of  the  system,  one  of  the  constraints 
of  the  problem  is  that  both  c  ,  the  number  of  repair  channels, 
and  y  ,  the  number  of  spares,  are  integer.  Therefore,  in  Figure 
I1I-1  we  are  interested  only  in  the  lattice  of  points  associated 
with  the  integer  pairs  (c,y)  .  The  solid  step-like  line  is  the 
boundary  of  the  region  for  which  the  pair  (c,y)  satisfies  the 
availability  constraint.  As  will  be  proven,  all  integer  pairs 
to  the  right  and  above  this  line  satisfy  the  availability  con¬ 
straint.  To  the  left  or  below  this  line,  the  constraint  is  not 
satisfied.  For  a  given  y  ,  the  minimum  c  for  which  the  avail¬ 
ability  constraint  is  satisfied  is  given  by  a  point  on  this  line. 

If  a  c  is  less  than  the  leftmost  value  of  the  line,  the  avail¬ 
ability  constraint  cannot  be  satisfied  for  that  c  ,  and  for  a  y 
less  than  the  lowest  value  of  the  line,  the  availability  constraint 
cannot  be  satisfied  for  that  y  . 

The  dashed  line  is  the  boundary  of  the  convex  hull  of  the 
region  satisfying  the  availability  constraint.  Each  of  the  ver¬ 
tices  of  this  convex  hull  must  be  at  one  of  the  (c,y)  pairs  for 
the  boundary  of  the  region  satisfying  the  availability  constraint. 
Not  all  of  the  (c,y)  pairs  that  lie  on  the  boundary  of  the  region 
satisfying  the  availability  constraint  will  lie  on  its  convex 
hull. 


The  solid  diagonal  line  c=y  is  important  in  this  enu¬ 
meration  technique  in  that  the  equations  specifying  the  state  of 
the  system  depend  on  which  side  of  the  line  the  state  falls. 
Equations  (II-l),  a,  b,  and  a  hold  when  c <  y  (to  the  left  of  the 
line),  and  Equations  (II-l),  d ,  e,  and  f  hold  when  c>y  (to  the 
right  of  the  line). 
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The  enumeration  procedure  starts  by  finding  the  smallest 
c  for  which  the  availability  constraint  is  satisfied,  and  deter¬ 
mines  the  minimum  y  associated  with  that  c  for  which  the  avail¬ 
ability  constraint  is  satisfied  (point  a).  Then,  as  will  be  shown, 
since  all  points  above  and  to  the  right  of  that  point  (a)  satisfy 
the  availability  constraint,  we  search  for  the  next  point  to  the 
right  and  below  for  which  the  availability  constraint  is  satis¬ 
fied.  The  procedure  goes  from  point  to  point  in  this  manner.  If 
we  are  searching  for  a  point  to  the  left  of  the  c=y  line.  Equa¬ 
tions  (II-l),  a,  b,  and  a  are  used.  As  soon  as  c>y  ,  Equations 
(II-l),  d,  e,  and  f  are  used.  The  procedure  continues  by  incre¬ 
menting  c  or  decrementing  y  until  the  availability  constraint 
cannot  be  satisfied  either  for  y=l  or  for  c  =  y+M  . 

It  is  proven  in  this  section  that  the  hull  can  be  generated 
in  this  manner  in  a  finite  number  of  steps.  The  outline  for  the 
algorithm  used  in  the  computer  program  is  also  presented. 

The  expression  given  in  Appendix  B,  (B-l),  for  the  avail¬ 
ability  constraint  will  be  used  as  we  examine  the  nature  of  the 
region  satisfying  the  availability  constraint  and  generate  expres¬ 
sions  used  in  the  computer  algorithm. 

By  the  nature  of  the  problem,  one  intuitively  feels  that 
if  the  availability  constraint  is  satisfied  for  a  (c,y)  pair, 
it  should  also  be  satisfied  for  (c',y’)  when  c'^c  and  y'^y  . 
This  shall  now  be  proven. 

Theorem  2:  If  the  pair  (c,y)  satisfies  the  availability  con¬ 
straint,  then  any  pair  (c',y')  ,  where  c'  >c  and  y'^y  , 
satisfies  the  availability  constraint. 

Proof:  Part  I}  a  <y  . 

The  state  probabilities  are  given  by  Equation  (II-l), 
a,  b ,  and  a,  which  are  repeated  here. 
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P 


n 


p0 


"5 

(0£nj<c)  a 


(c  <  n  <y) 


b 


(II-l) 


MyM!  n\n 

(M-n+y)!  cn-C  c!  P° 


(y  ^  n ^  y+M) 


a  . 


Using  Equation  (B-l)  from  Appendix  B,  the  availability  constraint  for 
c^y  can  be  written: 


c-1  n  y-1  n  ,  v  M-l  . ..  ,  /  \j 

y  Sl  +  y  — e —  >  /__§_ \  Dy  y  H-IHnI _ I _ /p  V 

n!  n-c  ,  -  la-l/P  M(M-j) !  j+y-c  ,  ImJ  » 

n=0  n=c  c  c!  '  '  j=0  J  c  c!\/ 


MA 

where  p  =  ,  and  p  has  been  eliminated  from  all  terms.  This 


p  ’  —  *0 

can  be  rewritten  as 

c-1  n  y-1  n 

y  £-r  +  y  — £ — 

n=0  n*  n=c  cn"C  c! 


/  v  M-l 

(*H  5, 


(M-l) ! 


j 


j=o  (M_1"j)!  ci+y-c  c!  W 

(in-3) 


Part  1(a) . 


It  will  now  be  proven  that  if  (III-3)  is  satisfied  for 
a  (c,y)  pair  then  it  will  be  satisfied  for  a  (c',y)  pair,  where 
c'  >c  and  c'  ^y  .  (It  will  subsequently  be  shown  that  the  avail¬ 
ability  constraint  remains  satisfied  for  c*  >y  .)  For  c'  =  c+1  , 
the  availability  constraint,  as  expressed  by  (III-3),  is 


y-1 

L 

n=c+l 


(c+l)n-c-1  (c+1)! 


> 


/_a_\oy  y1  --(M-l) ! _ 1 _ 

\1_a/  j=0  (M-1-!)1  (c+l)^-0'1  (c+1)! 


1 


or 
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c-L  n 

T  £t  + 


c 

£_  . 
c! 


y-i 


c: 


*—•  ,  , , .  n-c  , 

n=*c  (c+1)  cl 


(M-l) ! _ 1 

(M-i-j)!  (c+1)j+X'C  c| 


Multiplying  by  [(c+l)/c]y  c  ( 

ky-c  c-l  n  y-1 


/  .Ay-c  c-l  n  y-1  n  /  x1\y-n 

£  Sr  ♦  E  7^  (*)  i 

/_jJL y  r1  -Qtill _ 1 _ (py  uy 

\l-*r  h  (M-l-j)!  cj4y-c  c!  (mJ  (c+lj  • 


(111-4) 


Since  for  yj^c  ,  [(c+l)/c]y  C>1  ,  and  for  n<y-l  ,  [(c+l)/c]y  °> 

1  ,  a  comparison  of  the  left-hand  terms  of  (III-4)  to  the  left-hand 
terms  of  (III-3)  shows  that  both  of  the  terms  on  the  left-hand  side 
of  (111-4)  are  greater  than  the  equivalent  terms  in  (IH-3).  Also, 

since  for  j  >;  0  ,  (c/(c+l)]^<l  ,  the  right-hand  side  of  (II1-4)  is 
less  than  the  right-hand  side  of  (III-3).  Therefore,  if  the  relation 
(III-3)  is  satisfied,  then  the  relation  (III-4)  must  be  satisfied. 

So  if  the  availability  constraint  as  given  in  (III-3)  is  satisfied 
for  a  (c,y)  pair,  then  it  must  also  be  satisfied  for  (c+l,y)  , 
and  consequently  for  (c',y)  when  c'^c  . 


Part  1(b). 

Still  assuming  that  c^y  ,  it  will  now  be  proven  that 
if  (III-3)  is  satisfied  for  a  (c,y)  pair  then  it  will  be  satis¬ 
fied  for  (c,y')  when  yT  >y  . 


Let  y'  =  y+1  .  The  availability  constraint  can  then 

be  written 


c-l  „n  y  n  / 

— \py+1  r1 

(M-j)  M!  1  /p 

v 

n!  ^  n-c  *  =  (] 

n-0  n=c  c  c!  \ 

j?o 

(M-j)!  M  c j+y+l-c  c,  \M 

/  • 
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Rewriting  the  right-hand  side, 


y*1  pV  A  Pn  >  /  a  \  p  y  y^  (M-l) !  1  /p\j 

n=0  n!  n=c  cn"C  c!  =  V1’3/  C  h  (M_1"j)!  cJ4y_C  c!  W 


The  left-hand  side  can  be  written  as 


c— 1  n 


1  + 


n=l 


Regrouping  terms,  this  becomes 


y  2-  +  £_  +  y  — £ 

A  n!  c!  A,  n- 


n=c+l  c 


C  P° 


n=l 


Factoring  out  p/c  ,  we  get 


1  +  E  E  -S? 


n=c+l  c  c! 


c  n-1  y  n-1 

£+Ap__^+  y  — P - 

p  n!  n-l-c 


n-1 


n=c+l  c 


Redefining  the  indices  gives 


c-1  m  y-1 

£  +  y  £_  _£_  +  y 

o  m!  m+t 


m 


m=0 


m=c  c 


m-c 


cl 


therefore,  the  inequality  becomes 


-  ,  c-1  n 

£  £  +  y  2-^  +  y  — e 

c  )  p  A  n !  n+1  £->  n- 


y-i 


(*)« 


n=0 

M-l 

py  y  (m-1)  ! 


n=c  c 


C! 


j=o  c! 


Eliminating  the  p/c  factor  from  both  sides,  we  are  left  with  the 
inequality 

c  y"1  P^  c  A1  Pn  >  /  a  \  y  A*1  (M-l) ! _ 1  /pV 

P  n=0  n!  n+1  n=c  cn_C  c!  =  V1"3/  j=0  cj4y"C  ci  W 

(III- 5) 

Since  for  n  <  c-1  ,  {c/(n+l)]^l  ,  the  left-hand  side  of  (III-5)  is 
greater  than  the  left-hand  side  of  (III- 3).  Moreover,  the  right-hand 
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sides  of  (III-3)  and  (III-5)  are  identical.  Therefore,  if  (III-3) 
is  satisfied  for  a  (c,y)  pair  then  (III-5)  must  be  satisfied  for 
(c,y')  when  y'  =  y+1  ,  and  consequently  for  any  pair  (c,y*) 
when  y '  y  . 

In  Parts  1(a)  and  1(b)  of  this  proof  we  have  shown  that 
when  c<y  ,  if  a  pair  (c,y)  satisfies  the  availability  constraint 
then  (c',y)  and  (c,y')  satisfy  the  availability  constraint, 
respectively  when  c'  >c  and  c*  <y  ,  and  when  y'^y  .  This  im¬ 
plies  that  any  pair  (c',y')  ,  c’j>c  and  y'^y  ,  will  also  satisfy 
the  availability  constraint  when  c'^y  .  This  concludes  Part  I  of 
this  proof. 


Part  II,  a  . 

For  c  >y  the  state  probabilities  are  given  by  Equation 
(II-l),  d,  e,  and  /,  which  is  repeated  here. 


>>0  •  (Olni,)  d 

<«•*«>  • 

- -*1- -  •  (ff  P0  ■  (C£n<y«>  f 

(M-n-ty) !  c  c!  '  • 


pn  =  1  (M-n+y) 


(II-l) 


(M-n+y) 


Using  Equation  (B-l)  the  availability  constraint  for  c^y  can  be 


written  as: 


y-1  n  /v  c.-y- 

E  fr  i  A y  E 

n=0  '  !  j=0 


c  V'~1  (»~i)  H!  (pV 

A  (M-j)!  M( y+j ) !  \M/ 


+  y1  (M-j)  M!  /p\j 

i=c-y  (M-j)!  Mc"^  c  c! 


where,  again,  p  =  —  and  p^  has  been  eliminated  from  all  terms. 
This  can  be  rewritten  as: 


t 


V-1  >  l  a  V-,y  T  (M-l) !  /p\J 

n!  =  \l-a/P  j^0  (M-l-j ) !  (y+j) !  \M/ 

,  y,1  (M-l)  !  /p\j  ■ 

c!  j=r-y  <M-1'J>!  VM/ 


-(i+y) 


(III-6) 


Part  11(a). 


It  will  now  be  proven  that  if  (III-6)  is  satisfied  for 
a  (c,y)  pair  then  it  will  be  satisfied  for  a  (c',y)  pair  where 
c'_^c  .  For  c'  =  c+1  ,  (III-6)  becomes 


%  £  i  £ 

n=0  '  '  j=0 


(M-l)!  /p\J 

(M-l-j)!  (y+j)!  \yi) 

(IH-7) 


+  l£±l)__  V  (M-l)!  /pV  ,  ,n 
<c+l>!  j=c?y+l  W  (  } 


By  changing  the  limits  on  the  right-hand  summation  and 
eliminating  a  (c+1)  from  the  numerator  and  denominator  of  the  co¬ 
efficient  of  the  last  summation,  we  get 


yfV  >  /_&V  IT1 _ fltui _ M3 

ba  n!  -  M  I  jTo  (y+j)!  w 

+  (M-l) !  fp\C~y  (M-l)! _  /p\C' 

(M-l-c+y) !  c!  ^M/  (M-l-c+y) !  c!  \M/ 

+  <s*2)c-  V  (M-1)J_  /p\j  (c+1)-(j+y)l 

c!  .^_y  (M-l-j)!  [v)  W)  j  ’ 


y  £_  >  /_A_\py  |  y  _ (m-1)*  ..  /pV 

n=0  n!  ’  \1_a/  j-0  (y+J)»  W 

(II 

+  ^  j Iy  ($)’  <~»-<w  • 


(III-8) 


.  (c+l)C  y.1  (M-l) 

r !  CM-l-r 


c!  ,*-•  (M-l- 

j=c-y 
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Comparing  (III-8)  and  (III-6)  we  see  that  all  terms  are  identical 
except  for  the  last  summation  on  the  right.  Each  term  in  the  last 

summation  of  (III-8)  differs  by  a  factor  of  [c/(c+l) C  from 
the  comparable  term  in  (III-6). 

Since  in  the  last  summation  j  >  c-y  and  c/(c+l)  <  1  , 
it  is  obvious  that  the  last  summation  in  (III-8)  is  less  than  the 
last  summation  in  (III-6).  Therefore,  if  the  inequality  (III-6) 
is  satisfied,  then  the  inequality  ( I II— 7)  is  satisfied.  So  for 
c^y  ,  if  the  availability  constraint  is  satisfied  for  (c,y)  then 
the  constraint  is  satisfied  for  (c',y)  when  c'  =  c+1  and, 
therefore,  when  c'^c  . 


Part  II (b). 

It  will  now  be  proven  that  if  the  availability  con¬ 
straint  is  satisfied  for  a  pair  (c,y)  then  it  will  be  satisfied 
for  a  pair  (c,y')  when  y'^y  .  This  was  shown  to  be  true  when 
c  =  y  in  Part  1(b).  Now  we  shall  assume  that  y<c  and  will  show 
that  (III-6)  remains  satisfied  provided  that  y’j^c  . 


Let  y'  =  y+1  .  The  availability  constraint  can  then 


be  written 


y'-i 

„  \  „  f c— y 

E 

St  i  (i 

e 

n=0 

/  J-0 

(M-j) 


M! 


v-1  (M-1) _ M!  /p\j 

j4V  (M-j)!  M(c^'-C)  c!  W 


(HI-9) 


Substituting  y'  =  y+1  and  manipulating,  this  can  be  written 


£  *  *  (-)' 


,y+1 


^  c-y-2 

I  j50  "(M-1- j  )  !  (y+l+j  )  ! 


(M-1)! 


cl  y.1  (M-1)!  /pV  -( j+y+1) 

c!  j=c-y-l  (M"1-j)!  W 
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By  pulling  terms  out  of  the  summations  and  dividing  both  sides  by 
p  ,  the  inequality  can  be  rewritten  as 


y-1  m 


p  '  (m+1) 


.  M  y  1  _ (M~l)  ! _  Mj 

!  -  \l-a/p  j  (M-l-j ) !  (y+l+j) ! 

_  (M-l)!  M 

(M+y-c)!  c! 

.  c!  (M-l) !  /pV 
c!  (M+y-c)!  \MJ 

.  c!  v1  !  M 

C!  4^-y  (M-l-j)-'  \Hj 


kc-y-l 


,c-y-l 


i  c-(i+y) 


Multiplying  by  y  we  get 

y-i 


P 


E'  £_  _x_ 

.  n !  n+-l 
n=0 


,(4 


(M-l) ! 


c-y-1 

j5)  (M-i-j)-'  (y+j>!  y+l+j 


f  (")J 


(111-10) 


.  cC  v*1  (M-l)!  / pV  -(y+j)  y. 

c!  i=c-y  (M~1“i>!  W 


A  term  by  term  comparison  of  (111-10)  and  (II1-6)  shows 
that  if  the  inequality  (III-6)  is  satisfied  then  (III-10)  is  also 
satisfied.  On  the  left-hand  side  of  (III-10)  there  is  an  additional 
term  which  is  greater  than  zero  and  since,  in  the  left-hand  summa¬ 
tion,  n  <  y-1  ,  each  term  of  the  left-hand  summation  of  (111-10) 
is  greater  than  the  corresponding  term  in  the  left-hand  summation  of 
(I1I-6).  In  the  first  summation  on  the  right-hand  side,  since 
j  >  0  ,  each  term  of  (III-10)  is  smaller  than  its  corresponding  term 
in  (III-6),  and  in  tne  second  summation  on  the  right-hand  side,  since 
y  <  c  ,  each  term  of  (III-10)  is  also  smaller  than  or  equal  to  the 
corresponding  term  of  (III-6). 


Therefore,  if  for  a  pair  (c,y)  (1II-6)  is  satisfied, 

then  for  a  pair  (c,y+-l)  the  inequality  is  also  satisfied  and  the 
same  is  consequently  true  for  (c,y')  when  y'  >y  .  In  Parts  11(a) 


t 


A3 


and  11(b)  of  this  proof  we  have  shown  that,  for  yf-c  ,  if  a  pair 
(c,y)  satisfies  the  availability  constraint  then  a  pair  (c',y') 
also  satisfies  the  constraint  when  c'  >c  and  y'^y  . 

By  Parts  I  and  II  we  have  proven  that  for  y > c  or 
y<c  if  the  availability  constraint  is  satisfied  for  a  pair  (c,y) 
it  is  also  satisfied  for  (c',y’)  when  c'^c  and  y'^y  .  // 


Referring  back  to  Figure  III-l,  the  geometric  interpretation 
is  that  any  point  above  and  to  the  right  of  a  boundary  point  satis¬ 
fies  the  availability  constraint.  This  property  will  be  used  later 
in  the  generation  of  the  feasible  region. 


For  the  algorithmic  development  there  is  a  more  convenient 
form  for  (II1-3).  Note  that  the  second  summation  is  a  geometric 
sum.  Therefore,  for  p  ^  c  ,  (III-3)  can  be  written 


-1 

£ 

=0 


(*)  s  r  s  (43 


(M-l)  ! 
(M-l-j)!  * 


The  airantage  of  this  form  is  that  the  dependence  on  y  has  been 
remov™  from  inside  the  summation  and  from  its  indices.  It  remains 
only  as  an  exponent  of  the  term  p/c  . 


For  convenience,  the  following  notation  is  introduced: 

c-1  n  • 


si-£,f 

n=0 


V1  /_P_V  (M-l)  ! 

jT0  W  * 


It  is  obvious  that  and  are  both  nonnegative  and  finite. 

Also,  since  a  is  the  availability  (0<a<l)  ,  A  is  nonnegative 
and  finite. 


I 
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Using  this  notation,  the  availability  constraint  for  y  .>  c 
and  p?^c  can  be  written 


c  ulf'  c  /v 

+  -  ^  >  A  s2(t) 


1  c!  1_p 
c 


y-c 


or 


si  _J_ 


S1  +  cT  x_£  =  I  AS2  + 

c  \  c 


*)«  r  • 


Multiplying  by  c!/p  yields 


c ! 


S,  +  — 

PC  1  i-£  = 


>  fAS2  + 


y-c 


(III-ll) 


This  form  of  the  availability  constraint  is  very  useful  in  deter¬ 
mining  the  convex  hull  of  the  region  satisfying  the  availability 
constraint. 


For  c  >  p  ,  (111-11)  can  be  written 
ky-c 


<  I  ~  S1  + 

p  L  1-* 


P|/(AS2  +  ,  P 


1  -- 

C/ 


Taking  the  natural  log  of  both  sides, 


$ 


C ! 


(y-c)  £n[— J  <  £n  [  — 1  S  + 

pc  1  i  -  ^ 


Pi  -  £n  AS2+-p 


C/  \  C  y 

Since  (p/c)  <  1  ,  £n(p/c)  <  0  ;  therefore,  dividing  by  £n(p/c) 

gives 


(y- 


c)  ’  [*"  S1  +  7T|)  *  *"  (“2  +  TTi)]  ' tn© 


or 


+  L4l,  +  — i-' 

l  1/  1  i-f, 


-  In  I  AS2  + 


1  -■ 


c/J 


/  Zn(-  .  ( III-12) 


1 

I 

1 


* 


For  a  specified  c  ,  and  S2  can  be  evaluated.  The 

boundary  of  the  region  satisfying  the  availability  constraint  for 
a  given  value  of  c  is  at  the  minimum  integer  value  of  y  that 
satisfies  the  availability  constraint.  Let 

y.  ,  c  +  +  S1  +  TTp)  -  *»  • 

Then  for  a  given  c  ,  a  boundary  point  on  the  region  satisfying  the 
availability  constraint  is  given  by  the  pair  (c,Ty*l)  ,  where  fal 
is  the  smallest  integer  great  A:  than  or  equal  to  a  .  Remember  the 
relationship  (TII-12)  is  only  valid  for  y ^c  and  c>p  . 

For  y*  to  be  greater  than  or  equal  to  c  ,  since  £.n(p/c)  < 
0  ,  the  following  relation  must  hold: 

£  £"PS’+*)  ’ 


—  sn  . 

c  1 


Therefore,  if  (p  /c!)AS2  >  ,  we  know  that  for  a  specified  c 

the  point  (c, Ty*l)  is  a  boundary  point  of  the  feasible  region. 

If,  however,  >_  (pC / c!) AS ^  ,  we  know  from  (III-ll)  that  the  avail 

ability  constraint  is  satisfied  for  y=c  and  therefore,  either 
(c,y)  is  a  boundary  point  when  c  =  y  ,  or  (c,y’)  ,  where  y'  <c  , 
is  a  boundary  point.  If  (c,y')  when  y'  <  c  is  the  boundary 
point,  another  method,  as  described  later,  is  needed  to  determine 
that  boundary  point. 


For  c  <  p  and  y_>c  ,  we  can  derive  another  relationship 
from  (III-ll).  For  c<p  ,  l/[l-(p/c)]  <  0  and  therefore, 

AS^  +  1/[1-  (p/c)]  may  or  may  not  be  greater  than  zero  and 

(c!/pc)S^  +  1/ [ 1  —  (p/c)]  may  or  may  not  be  greater  than  zero.  How¬ 


ever,  we  can  prove  that 
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—  S  +  — - —  <  0 

pC  1  l-£  = 

c 


Theorem  3:  For  c<p  and  y>c  , 


—  S  +  — <  0  . 

pc  1  1-fi  = 


(III-13) 


Proof:  The  definition  of  is 


Si  =  n=0 


c-1  n 

iESt-  *>° 


We  can  write  (III-13)  as 


S.  < 

1  =  c! 


PC  1 


Since  p>c  ,  the  right-hand  side  is  obviously  positive.  For  c  =  0  , 
is  defined  to  be  equal  to  0,  so  the  inequality  if  satisfied. 

For  c > 1  ,  the  right-hand  side  can  be  expanded  as  a  geometric  pro¬ 
gression: 


c  °°  /  \  ra  c-1  /  \  m 


Using  the  definition  of  and  splitting  the  summation 


on  the  right,  we  get 


c-1  n  c-1  m 

y  £_  <  y  — c_ 

“  n!  --  (c-1 

n=0  ra=0 


c-1  m  .  c-1  °°  ,  »m 

y  —5 _  c-l-m  _p -  y  (c\ 


By  changing  the  index  of  the  first  sum  on  the  right  (n=  c-l-m)  and 
bringing  it  over  to  the  left,  the  relation  can  be  written 


m 
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c-l  (  ,  c-l-n  c-l  ®  /  \n 


(III-14) 


We  will  now  show  that  the  left-hand  side  of  this  in- 
equality  is  negative.  We  must  show 

,  c-l-n 


n!  (c-l) ! 
is  negative;  i.e.,  we  must  show  that 

c-l-n 


1  ft  n!  c-l-n\ 
n!  "  (c-l) !  c  ) 


n!  c 


(c-l) ! 


>  1  , 


which  is  equivalent  to 


(c-l) (c-2)  ...  (n+1) 


>  1  , 


which  is  obviously  true.  Thus  the  left-hand  side  of  (^LII-14)  is 
negative.  Th#  sum  on  the  right-hand  side  is  obviously  positive,  so 
that  the  inequality  holds.  // 


Therefore,  since  (c!/pC)S,  +  1/[1-  (p/c)]  <  0  ,  if  the  avail- 

i  = 

ability  constraint  is  to  be  satisfied  for  a  real  y  ,  ASj  + 

1/[1-  (p/c)]  also  must  be  less  than  or  equal  to  zero.  For  AS^  + 

1/ [ 1  —  (p/c)]  <  0  ,  (III-ll)  can  be  written 

(S'** ri|) ri)- 

Sin^e  p  >  c  ,  taking  the  natural  log  of  both  sides  and  dividing  by 
£n(p/c)  ,  we  get 

(y_c)  >^[ps1+_y  /  j“2+^|J]  //"(f). 


or 
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ft 


>  c  +  Jin 


This  is  an  equivalent  form  of  (III-12);  therefore,  again  the  bound¬ 
ary  of  the  region  satisfying  the  availability  constraint  for  a  given 
value  of  c  is  at  the  minimum  integer  value  of  y  that  satisfies 
the  availability  constraint.  Let 


y* 


c  +  Jin 


Then  for  a  given  c  ,  a  boundary  point  of  the  region  satisfying 
the  availability  constraint  is  given  by  the  pair  (c,fy*l)  .  But 
y*  must  be  greater  than  or  equal  to  c  for  (III-ll)  and  the  subse¬ 
quent  analysis,  to  be  true.  Therefore,  since  £n(p/c)  >  0  , 


Jin 


must  be  greater  than  or  equal  to  zero.  This  implies  that 

c  i  1  1 

—  S.  +  -  >  AS.  +  , 

pc  1  i-£  “  2  i-£ 


or 


ci 


c  S1  *  As2  ’ 


thus. 


si  *  cT  as2 


But  this  is  the  availability  constraint  for  y *  c  .  Therefore,  for 
c<p  ,  either  the  availability  constraint  is  satisfied  for  y<c 
or  it  is  not  satisfied  at  all.  This  property  is  useful  in  develop¬ 
ing  the  algorithm  for  generating  the  boundary  of  the  region  satisfy¬ 
ing  the  availability  constraint,  which  is  presented  below. 


Finally,  for  y>c  ,  there  is  one  more  case  of  interest,  and 
that  is  when  p  =  c  .  For  this  case  expression  (III-ll)  is  not  valid. 
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Instead,  using  (III-3),  we  obtain 


c-1 

£ 

n=0 


> 


a 

1-a 


p!  Zlf.  v1  o1-1)  1  LpY 

c!  cy-c  £0  0*-J)!  VMc/  ‘ 


V-c  v— c 

Since  p  =  c  ,  (p  /cJ  )  »  1  and  the  second  summation  on  the  left 
equals  (pC/c!)(y-c)  •  Substituting  A  ,  ,  and  S^  for  the  ap¬ 

propriate  terms,  we  get 


—  S.  +  y-c  >  AS- 
cl  “  2 

P 

or 

y  >  c  +  AS-  -  —  S..  . 

2  pc  1 

Let 

y*  =  c  +  AS  -  —  S  . 

2  pc  1 


Then  for  a  given  c  ,  a  boundary  point  of  the  region  satisfying  the 
availability  constraint  is  given  by  the  pair  (c,ry*l)  .  As  for  the 
case  c  >  p  ,  y*  will  be  greater  than  or  equal  to  c  if 


which  is  equivalent  to  the  situation  in  which  the  availability  con¬ 
straint  S2  >.  (pC/c!)AS^  is  not  satisfied  for  y=c  .  Therefore, 

for  p=c  ,  either  the  availability  constraint  is  satisfied  for 
y=c  or  there  is  a  y'  >  c  such  that  (c,y')  is  the  boundary 
point  for  the  region  satisfying  the  availability  constraint  for  a 
specified  c  . 


For  y<c  ,  a  convenient  form  of  the  availability  constraint, 
similar  to  (III-ll)  for  y  >  c  ,  has  not  been  found. 

From  all  of  the  above,  it  is  now  possible  to  state  an  algo¬ 
rithm  that  will  identify  all  of  the  (c,y)  pairs  that  define  the 


aaBMHiBCiaMauBta 
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boundary  of  the  region  satisfying  the  availability  constraint  for 
a  given  time  period. 

Algorithm. 

1.  For  the  time  period  of  interest,  calculate 
p  =  M(X/u)  . 

2.  Start  with  c  =  fplt. 

3.  Calculate  AS2  ,  where  A  =  a/(l-a)  and  = 

Y'1  /fiV  (M-l)  ! 

JT0  VMc/  * 

c-l  n 

4.  Calculate  S  ,  where  S  =  —  . 

1  1  n=0  "* 

5.  If  (c!/pC)S1  <  AS2  ,  go  to  Step  12. 

6.  If  c>  f pi  ,  go  to  Step  18. 

7.  Decrement  c  by  1  . 

8.  Calculate  A  ,  S^  ,  Sj  • 

9.  If  (c ! / pC) S^  >  AS 2  ,  go  to  Step  7. 

10.  Increment  c  by  1  . 

11.  Go  to  Step  18. 

12.  If  c 4  P  ,  go  to  Step  15. 

13.  y  =  fc  +  AS2  -  (c!/pC)S11  . 

14.  Go  to  Step  16. 


tfal  is  the  smallest  integer  greater  than 


a  . 
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15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


25. 


The 

explained . 


1. 


2. 


(c,y)  is  a  boundary  point.  Increment  c  by  1  . 
Go  to  Step  3. 
y  =  c-1  . 

Test  if: 


n=0 


n 

( c-y-1 

fr  i  * 

1  T 

!  l-o 

(M-l 

c 

M-l 

+  rr . 

E 

3 

=c-y 

(M-l) ! 


If  the  inequality  is  satisfied,  go  to  Step  23. 

The  point  (c,y+l)  is  on  the  boundary.  Increment 
c  by  1  and  continue  the  algorithm. 

If  c  =  y+M+1  ,  the  algorithm  is  terminated  and 
all  boundary  points  have  been  generated. 

If  y  =  1  ,  terminate  the  algorithm.  All  boundary 
points  have  been  generated. 

Decrement  y  by  1  . 

Go  to  Step  19. 


purpose  of  each  step  of  the  algorithm  will  now  be 


Calculate  p  ,  a  parameter  of  the  problem. 

Initialize  c  equal  to  the  greatest  integer  less 
than  p  .  If  there  is  a  boundary  point  with  y  >  c  , 
then  c  >  p  . 
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3.  Calculate  AS£  •  a  term  needed  to  determine  if 

the  availability  constraint  is  satisfied  for 
y  >  c  . 

4.  Calculate  ,  another  term  needed  to  evaluate 
the  availability  constraint  for  y  >  c  . 

5.  Test  if  the  boundary  point  is  such  that  y  >  c  . 

6.  The  next  boundary  point  has  y <  c  .  Test  if  this 
is  the  first  boundary  point  generated. 

7.  This  is  the  first  boundary  point  generated. 

Search  for  the  minimum  c  such  that  the  avail¬ 
ability  constraint  is  satisfied. 

8.  Calculate  needed  terms. 

9.  If  the  availability  constraint  is  satisfied  for 
y<c  ,  then  it  will  be  satisfied  for  y=c  . 

Test  to  see  if  the  availability  constraint  is 
satisfied  for  y“c  .  If  it  is  still  satisfied, 
go  to  Step  7  to  continue  searching  for  the  first 
c  for  which  the  availability  constraint  is  not 
satisfied. 

10.  The  first  c  for  which  the  availability  con¬ 
straint  is  not  satisfied  has  been  determined. 

Set  c  to  the  last  c  for  which  the  availabil¬ 
ity  constraint  can  be  satisfied. 

11.  Go  to  Step  18  to  determine  the  boundary  when 
y  <  c  . 

12.  Test  if  c  equals  p  . 

13.  Calculate  the  boundary  y  when  c  =  p  . 

14.  The  boundary  point  is  identified. 

15.  Calculate  y  for  the  boundary  point  when  c  > p 
and  y  >  c  . 


16. )  The  boundary  point  is  determined  for  y^c  . 

17. )  Proceed  to  calculate  the  next  boundary  point. 

18.  Boundary  point  calculation  for  y£c  . 

19.  Determine  whether  the  availability  constraint  is 
satisfied. 

20.  Test  if  availability  constraint  is  satisfied.  If 
it  is  satisfied,  continue  to  decrement  y  for 
the  specified  c  to  determine  the  minimum  y  . 

21.  The  availability  constraint  is  no  longer  satis¬ 
fied  for  the  particular  c  .  Therefore,  the  pre¬ 
vious  (y,y+l)  was  the  minimum  y  for  the 
specified  c  for  which  the  availability  con¬ 
straint  was  satisfied.  Proceed  to  the  next  c  . 

22.  If  c  is  greater  than  y+M  we  know  that  there 
are  no  more  boundary  points.  This  is  because  of 
the  physical  reality  as  well  as  the  mathematical 
formulation.  There  will  never  be  more  repair 
channels  than  there  are  machines  in  operation 
plus  spares. 

23.  We  will  never  have  negative  spares,  so  if  y  =  0 
we  will  not  have  any  additional  boundary  points. 

24.  Continue  the  search  for  the  minimum  y  that 

25. )  satisfies  the  availability  constraint  for  a 

specified  c  . 

The  algorithm  is  best  illustrated  by  Figures  III-2  and  1II-3. 
We  start  by  choosing  c  *  fpl  .  If  the  availability  constraint  is 
satisfied  for  some  y  less  than  or  equal  to  c  we  are  on  Figure 
HI-3;  otherwise  we  are  using  Figure  III-2.  For  Figure  III-2  we 
proceed  from  c  =  fpl  and  increment  c  by  1  ,  finding  the  minimum 
y  that  satisfies  the  availability  constraint  for  each  c  .  Eventu¬ 
ally  the  boundary  crosses  the  line  y=c  .  To  the  right  of  this 
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line,  using  Steps  18-25,  y  is  decremented  by  1  and  the  minimum 
c  that  satisfies  the  availability  constraints  for  the  specified  y 
is  determined  by  incrementing.  This  procedure  continues  until  ei¬ 
ther  y=l  or  the  availability  constraint  is  not  satisfied  for 
c  =  y+M  .  For  Figure  III-3  there  is  no  boundary  point  for  y  >  c  . 

We  proceed  from  c  =  Tpl  to  the  minimum  c  for  which  there  is  a 
y  such  that  the  availability  constraint  is  satisfied.  Once  this 
c  is  determined  we  proceed  in  the  same  manner  as  above  for  the 
cases  when  y  <  c  . 

This  procedure,  therefore,  is  one  in  which  the  minimum  c 
for  which  the  availability  constraint  is  satisfied  is  determined, 
and  by  incrementing  c  until  either  c  =  y+M  or  y  = 1  ,  the  mini¬ 
mum  y  associated  with  each  c  for  which  the  availability  con¬ 
straint  is  satisfied,  is  determined.  The  algorithm  starts  from  the 
upper  left  of  the  region  and  proceeds  until  the  boundary  crosses 
either  the  abscissa  or  the  line  c  =  y+M  .  The  initial  c  either 
equals  the  greatest  integer  less  than  p  or  lies  below  the  line 
y=c  .  Therefore,  if  p  is  finite,  the  initial  c  is  finite. 

If  AS2  is  finite,  then  the  initial  y  is  finite.  Therefore,  the 

boundary  of  the  region  satisfying  the  availability  constraint  can  be 
described  by  a  finite  number  of  (c,y)  pairs. 

Theorem  4:  If  p  and  AS2  are  finite,  the  boundary  of  the  region 

satisfying  the  availability  constraint  can  be  enumerated  by  a  finite 
number  of  (c,y)  pairs. 

Proof:  For  a  finite  p  ,  the  minimum  c  ,  say  c*  ,  for  which 

the  availability  constraint  is  satisfied  is  finite,  c*  <  p 
If  c*  “  p  and  AS2  is  finite,  the  value  of  y  ,  say  y*  , 

corresponding  to  the  boundary  is  given  by  (III-12)  and  is  finite. 

For  c*  <  p  we  have  shown  that  y  must  be  less  than  c*  .  There¬ 
fore,  the  y*  associated  with  c*  ,  the  minimum  c  for  which  the 
availability  constraint  is  satisfied,  is  finite. 
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We  have  shown  that  when  a  pair  (c,y)  satisfies  the 

availability  constraint,  then  the  pair  (c',y')  ,  where  c'  >c  and 

y'  >y  ,  also  satisfies  the  availability  constraint.  Therefore,  if  a 

4*  4*  4"  41 

pair  (c  ,y  )  is  on  the  boundary,  then  either  c  <  c  or  y  <  y  . 

But  <_*  is  the  minimum  c  for  which  the  availability  constraint  is 

satisfied.  Therefore,  for  any  other  boundary  point  satisfying  the 

availability  constraint,  it  must  be  true  that  c+ > c*  .  So,  for  all 

other  boundary  points  y+ < y*  . 

Then  we  have  y* > y  > 0  and  y  is  integer.  Since  y* 
is  finite  for  finite,  there  are  a  finite  number  of  points 

needed  to  enumerate  the  boundary  of  the  region  satisfying  the  avail¬ 
ability  constraint.  // 

The  use  of  this  algorithm  provides  a  method  of  enumerating 
the  boundary  of  the  region  satisfying  the  availability  constraint 
for  each  time  period.  We  know  that  for  a  given  c  (or  a  given  y  ) 
the  minimum  y  (c)  for  which  the  availability  constraint  is  satis¬ 
fied  is  on  the  boundary.  If  the  integer  constraint  and  the  avail¬ 
ability  constraints  were  the  only  constraints  of  the  problem  with 
a  linear  objective  function,  the  solution  would  be  on  this  boundary 
in  each  time  period.  There  are,  however,  two  additional  types  of 

constraints,  c.  >  c.  ,  and  y.  >  y.  ,  .  The  enumeration  scheme 

l  =  i-l  i  =  l-l 

must  account  for  these  constraints  as  well.  To  handle  these  con¬ 
straints  a  construct  called  a  "corner"  will  be  defined.  This  is 
described  next. 


III. 5  Corners 


In  this  section  we  will  define  a  quantity  which  we  term  a 
"corner"  and  describe  its  properties.  Through  the  use  of  corners, 
the  enumeration  of  the  overall  feasible  region  is  achieved. 


Definition:  For  a  given  point  a  =  (a^.a^, . . - ,am>  the  corner 
generated  by  a  is 

A  =  (x^a)  , 

where  x  =  (x,  ,x„,...,x  )  . 

1  c  ra 

Given  points  a^,...,aP  which  generate  corners  A^,...,AP  , 

we  define  a  dominated  region  or  D  region  as  D  =  U  A1  , 

i 

Thus  D  contains  all  points  contained  in  at  least  one  of  the 
corners. 

Given  any  region  R  of  the  form  D  above,  there  is  a  unique 

minimum  set  of  points  G(R)  which  generates  R  ;  i.e.,  a1 £  G(R) 

if  and  only  if  a1 e  R  ,  but  a1  is  not  dominated  by  another  point 
in  R  . 


Theorem  5:  Let  A1  and  A"1  be  corners  generated  by  a1  and  a^  . 

i  "i  1c  k  k 

Then  A  H  AJ  defines  a  corner  A  generated  by  a  ,  where  a^  = 

maxCa1  ,  a~*)  ,  n=l,2  . 
n  n 


1c  X  i 

Proof:  By  definition  A  =  {x^a  and  x_>aJ}  .  Therefore, 

k  x  *i 

if  xeA  ,  we  must  have  x  >  max{a  ,  aJ }  .  In  particular,  the 

n  =  n  n  r 

k  k  x  i  k 

point  a  defined  by  a  =  max(a  ,  aJ)  is  in  A  .  Clearly  no 

n  n  n 

k  k  k  k 

other  point  in  A  dominates  a  ,  so  that  a  generates  A 


Figure  III-4  illustrates  this  property  of  corners.  There  are 
two  separate  nresentations  given  in  Figure  III-4:  A  and  B.  In 

1  2 

both,  there  are  two  corners  generated  by  two  points:  a  and  a 

In  Figure  III-4A  we  see  that  a  D  region  is  generated  using  both 

12  12 
points  a  and  a  ,  and  the  intersection  AHA  is  generated 


by  a  new  corner  generated  by  a  •  In  Figure  III-4B,  the  D  region 

1  112 
is  the  same  as  the  corner  A  generated  by  a  ,  and  AHA  is 

2  2 
A  ,  the  same  as  the  corner  generated  by  a 

Another  property  of  D  regions,  which  will  prove  useful,  is 
the  manner  in  which  the  intersection  of  two  D  regions  can  be  gen¬ 
erated  by  the  generators  of  each  of  the  D  regions.  Let  = 

U  A*  and  D_  =  U  A^  .  Then 
iel  1  jeJ  ^ 


This  property  can  be  visualized  from  Figure  III-5.  The  dom- 

12  3 

inated  region  is  generated  by  the  set  of  comers  (A  ,A  ,A  ) 

12  3 

associated  with  the  points  (a  ,a  ,a  )  ,  and  the  dominated  region 

12  3 

D2  is  generated  by  the  set  of  corners  (B  ,B  ,B  )  associated  with 

12  3  ~ 

the  points  (b  ,b  ,b  )  .  The  new  region,  denoted  by  D  ,  is  gener- 

1  2  3  A  5 

ated  by  the  set  of  corners  (C  ,C  ,C  ,Cr)  generated  by  the 
1  2  3  A  5 

points  (c  ,c  ,c  ,c  ,c  )  ,  determined  via  Theorem  5.  Note  that  some 
1  2 

pairs  (e.g.,  a  and  b  )  produce  a  corner  which  is  contained  in 
some  other  corner  (e.g.,  ). 

The  utility  of  corners  arises  from  the  properties  of  the 
region  satisfying  the  availability  constraint,  as  discussed  in  the 
previous  section.  As  stated  in  Theorem  2,  if  (c',y')  satisfies 
the  availability  constraint  then  (c,y)  satisfies  the  constraint 
where  c>c'  and  y  >y'  •  Therefore,  if  the  point  (c’,y')  satis¬ 
fies  the  availability  constraint,  all  of  the  pairs  contained  in  the 
corner  generated  by  the  point  (c',y')  also  satisfy  the  availabil¬ 
ity  constraint. 
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For  any  given  time  period  i  ,  the  region  satisfying  the 
availability  constraint  can  be  described  by  a  D  region,  . 

The  boundary  of  has  the  property  that  for  (c^  ,  yj)  on  the 

boundary  of  then  either  (a)  for  y’  fixed,  cj  is  the 

smallest  integer  ^  for  which  (c^  ,  yp  satisfies  the  availabil¬ 
ity  constraint  for  time  period  i  ;  or  (b)  for  cj  fixed,  y^  is 
the  smallest  integer  y  for  which  (cj  ,  yp  satisfies  the  avail¬ 
ability  constraint.  If  both  properties  (a)  and  (b)  are  satisfied, 
then  (c^  ,  yp  generates  a  corner  composing  the  D  region  . 

The  enumeration  algorithm,  specified  in  the  previous  section,  gen¬ 
erates  the  unique  minimum  set  of  points  G(Dp  for  region  . 

We  have  defined  the  D  region,  ,  that  contains  all  of 

the  points  that  satisfy  the  availability  constraint  for  the  1th  time 
period.  We  have  also  specified  the  properties  of  the  boundary  of 
and  presented  a  method  of  enumerating  the  points  that  generate 

the  corners  composing  .  Besides  the  integer  constraint  and  the 

availability  constraint,  there  remain  the  monotonicity  constraints, 

i.e.,  c. >  c.  >  c.  ,  and  y. >  y.  >  y.  .  .  We  can  use  the 
i+1  =  i  =  i-1  l+l  =  i  —  i-1 

properties  of  D  regions  to  eliminate  points  that  do  not  satisfy 

these  constraints. 

Theorem  6:  If  the  pair  (c^,yp  satisfies  the  availability  con¬ 
straint  for  time  period  i  ,  then  for  (ci»yp  to  feasible  it 

must  also  satisfy  the  availability  constraint  for  time  period  j  , 
for  all  j  less  than  i  . 

Proof:  By  the  monotonicity  constraints,  c^  >  c^_^  .>  ...  >_ 

and  y±  >  y^  >  . . .  >  y^  . 
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From  Theorem  2  we  know  that  if  satisfies  the 

availability  constraint  in  time  period  j  ,  then  (c^  ,  yj)  , 
and  y^  >  y^  ,  must  also  satisfy  the  availability  constraint  for 
time  period  j  . 


Therefore,  since  c^  >  c^  ,  j  <  i  and  y^  >  y^  ,  j £ i  , 
then  the  pair  (c^,y^)  must  satisfy  the  availability  constraint 
for  the  time  period  i  for  all  j  <.  i  .  // 


Therefore,  if  the  point  associated  with  the  pair 
(c^,y^)  that  satisfies  the  availability  constraint  for  time  period 

i  is  feasible,  it  must  not  only  be  contained  in  ,  but  must  also 

be  contained  in  for  all  j  less  than  i  .  This  implies  that 

the  point  (c^.y^  is  contained  in  the  intersection  of  and 


i-1 


i.e.,  e  Di^Di  i 


Define  (there  is  no  time  period  previous  to  the 

first),  and  set  .  Then  for  a  point  (c^.y^  to  be 

feasible,  i.e.,  c^  and  y^  integer,  satisfies  the  avail¬ 

ability  constraint,  and  c^  and  y  satisfy  the  monotonicity  con¬ 
straints,  it  must  not  only  be  contained  in  but  it  must  also  be 

contained  in  .  Now  can  be  explicitly  determined  using  the 

property  of  the  intersection  of  two  D  regions  and  through  the  use 
of  Theorem  5. 

Let  us  further  examine  the  nature  of  the  boundary  of  . 


Then 


Theorem  7:  Let  (cVy)  be  on  the  boundary,  35^,  of  . 

(a)  for  y  fixed,  'c  is  the  smallest  integer  c 

for  which  (c.'y)  satisfies  the  jth  availability 
constraint  for  some  j  <  i  ;  or  , 

(b)  for  "c  fixed,  y  is  the  smallest  integer  y 

for  which  (c'.y)  satisfies  the  kth  availabil¬ 
ity  constraint  for  i  . 

Proof:  The  proof  is  by  induction. 

For  j  ■  2  ,  let  (c'.'y)  £  3D2  .  Then  either  (cVy)  e 
3D^  or  (c'.'y)  £  3D2  ,  since  I)2  *  D^nD2  .  Therefore,  (c’.'y)  sat¬ 
isfies  either  (a)  or  (b)  for  time  period  1  or  2,  respectively.  The 
theorem  is  therefore  true  for  j  =  2  (it  is  obviously  true  for 

j-1  >• 


Then  let  us  assume  it  is  true  for  i  =  j  .  We  must  show 
that  it  is  true  for  k  *  j+1  .  Let  (c'.'y)  e  •  Then  either 

(c'.y)  e  3°j+1  or  (~»y)  £  .  since  D^+1  -  D^nDj  *  If 

(cVy)  e  ^Dj+l  ’  then  (a^  or  ^  is  true  fey  c^e  definition  of  D^+1  • 

If  (~,*y)  £  3d"j  ,  the  theorem  holds  by  the  inductive 
assumption.  // 


Theorem  7  implies  that  a  point  (c^y^)  on  the  boundary  of 
Uj,  is  on  the  boundary  of  the  D  region  satisfying  the  availability 
constraint  ,  or  that  a  monotonicity  constraint  is  binding;  i.e.. 


ci-l  or  yi 


l-l  ' 


As  was  stated  previously,  a  linear  objective  function  will 
achieve  its  minimum  at  an  extreme  point  of  the  convex  set  of  all 
possible  solutions.  This  set  is  the  feasible  region  as  determined 
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by  the  constraints.  The  spares  provisioning  problem,  as  formulated, 
has  three  types  of  constraints:  the  integer  constraint,  the  avail¬ 
ability  constraint,  and  the  monotonicity  constraints.  The  extreme 
points  associated  with  the  integer  and  availability  constraints  are 
contained  in  the  set  of  points  generating  the  corners  of  the  D 
regions,  .  The  inclusion  of  the  backward  monotonicitv  constraint, 

i.e..  c,  >  c,  ,  and  y.  >  y,  ,  ,  is  accounted  for  in  the  generation 
1  =  i-1  i  —  i-1 

of  the  S’  .  There  remains  only  one  type  of  constraint  on  c^  and 

y^  that  has  not  been  accounted  for  in  the  generation  of  the  . 

These  are  the  forward  monotonicity  constraints,  i.e.,  c^  <  c^+^ 

and  y.  <  y. .  By  this  it  is  meant  that  if  c.  =  c. ,,  or  y.  = 

i  =  l+l  J  x  i+1  l 

y^+^  ,  there  is  no  assurance  that  the  associated  point  is  on  the 

boundary  of  .  This  last  constraint,  however,  can  be  accounted 

for  in  the  dynamic  programming  procedure,  as  discussed  in  the  next 
section. 


For  the  last  time  period  there  is  no  such  constraint.  There¬ 
fore,  the  extreme  points  for  the  last  time  period,  d  ,  must  come 
from  the  set  of  points  that  generate  the  corners  of  .  This  is  a 


finite  set.  For  all  other  time  periods,  j  <  £  ,  there  is  a  finite 
set  of  points  that  generate  the  regions.  These  points  contain 

the  set  of  extreme  points  that  satisfies  the  availability  constraint, 
integer  constraint,  and  backward  monotonicity  constraint,  i.e.. 


c  >  c 
i  -  j-i 


and  y  >  y 


j  -  'j-1  ' 


III. 6  Dynamic  Programming 


As  dynamic  programming  is  a  general  problem  solving  approach, 
there  does  not  exist  a  standard  mathematical  formulation.  However, 
there  are  certain  characteristics  associated  with  a  problem  if  dy¬ 
namic  programming  can  be  applied.  They  are: 


66 


1.  The  problem  can  be  divided  into  stages  with  a  pol¬ 
icy  decision  required  at  each  stage. 

2.  There  are  a  number  of  states  associated  with  each 
stage. 

3.  The  effect  of  the  policy  decision  at  each  stage  is 
to  transform  the  current  state  into  a  state  asso¬ 
ciated  with  the  next  stage. 

4.  Given  the  current  state,  an  optimal  policy  for  the 
remaining  stages  is  independent  of  the  policy 
adopted  in  the  previous  stages. 

5.  The  solution  procedure  begins  by  finding  the  opti¬ 
mal  policy  for  each  state  of  the  last  stage. 

6.  One  can  establish  a  recursive  relationship  that  re¬ 
lates  the  optimal  policy  for  each  state  at  the  cur¬ 
rent  stage  to  the  given  optimal  policy  for  each  state 
at  the  next  stage. 

7.  The  solution  procedure  moves  backwards,  stage  by 
stage,  until  the  .iptimal  policy  starting  with  the 
first  stage  is  found  [9]. 

At  first  glance  the  spares  provisioning  problem  appears  to 
satisfy  all  of  these  characteristics.  There  are  identifiable  stages 
associated  with  each  time  period.  The  states  are  associated  with  the 
choice  of  the  combination  of  the  number  of  spares  and  repair  channels 
at  each  time  period.  The  linear  objective  function  provides  a  rela¬ 
tionship  upon  which  to  determine  the  optimal  policy  at  each  stage. 

On  second  thought,  the  problem  becomes  formidable.  There  are 
two  major  obstacles  to  the  direct  application  of  dynamic  programming 
techniques  to  the  spares  provisioning  problem.  First,  the 
states  of  the  problem  cannot  be  specified  for  any  but  the  last  and 
next  to  last  stages  (time  periods  1  and  2),  because  the  X^’s  cannot 

be  determined  until  subsequent  stages  are  solved.  Secondly,  for 
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each  stage  there  is  potentially  an  infinite  number  of  states 
f(c,y)  pairs  that  satisfy  the  constraints]. 


I 


As  has  been  discussed,  A  for  i>2  is  dependent  on  the 


choice  of  c. 

1 


and 


where  j  <  i  .  In  a  standard  dynamic  pro¬ 


gramming  solution  the  problem  is  solved  working  backwards  in  time, 
i.e.,  the  first  stage  solved  is  associated  with  the  last  time  period. 
Therefore,  A^  ,  for  a  particular  time  period,  cannot  be  determined 


until  all 


pairs  are  specified  for  j  <  i 


Although  most 


dynamic  programming  solutions  do  proceed  backwards  in  time,  there  are 
techniques  for  proceeding  forward  [10],  These,  however,  require 
knowledge  of  the  initial  state  of  the  problem.  For  the  spares  pro¬ 
visioning  problem,  however,  there  is  potentially  an  infinite  number 
of  initial  states. 


This  first  problem,  the  indeterminateness  of  A^  ,  can  be 


resolved  by  using  the  branch  and  bound  procedure  presented  previous¬ 
ly.  Through  the  branch  and  bound  procedure,  subproblems  are  created 


in  which  the  X^'s  are  replaced  by  A^  .  The  A^  are  independent 
of  the  choice  of  (c^,y^)  for  time  periods.  Therefore,  it  is  pos¬ 
sible  to  determine  all  states  at  each  stage  of  the  problem,  even 
when  proceeding  backwards  in  time. 


The  second  problem,  the  number  of  potential  states,  can  be 
eliminated  by  the  enumeration  technique  and  the  use  of  Corners. 
There  is  an  infinite  number  of  choices  of  c  and  y  that  satisfy 
the  four  types  of  constraints:  c  and  y  are  integer;  the  avail¬ 
ability  constraint;  c±  <  ci+1  ,  yi  <  y±+1  ;  and  ci  >  ,  yt  > 

y^  ^  .  We  know,  however,  due  to  the  linearity  of  the  objective 

function,  that  the  solution  will  lie  on  the  extreme  points  of  the 
boundary  of  the  region  satisfying  these  constraints.  The  enumera¬ 
tion  technique  accounts  for  the  integer  constraint  and  determines 
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the  boundary  of  the  region  satisfying  the  availability  constraint. 

The  technique  associated  with  the  corners  modifies  the  boundary  due 
to  the  backward  monotonicity  constraints  (c^  >  c^^  ,  >  y  ^  . 

The  only  constraints  not  accounted  for  by  techniques  already  speci¬ 
fied  are  the  forward  monotonicity  constraints  (c^  <  c^+^  ,  y^  < 

y^+^)  which  are  accounted  for  directly  in  the  dynamic  programing 

procedure.  In  determining  the  set  of  extreme  points  on  the  boundary 
satisfying  the  first  three  of  the  four  constraint  types  for  each 
time  period,  we  end  up  with  a  finite,  and  usually  small,  set  of 
points. 

The  forward  monotonicity  constraints,  c^  £  Ci+i  »  £  ^i+1  * 

do  not  affect  the  use  of  dynamic  programming.  These  constraints  re¬ 
late  the  optimal  choice  of  c^  and  y^  in  a  particular  time  period 

to  the  selection  of  c  and  y  in  the  next  time  period.  Since 

1+1  1+1 

the  dynamic  programming  procedure  works  backward  in  time,  an  optimal 
choice  of  c^+^  and  Yi+1  is  made  before  the  optimal  c^  and  y^ 

are  determined.  Thus  once  the  optimal  policy  is  found  for  a  particu¬ 
lar  stage  (associated  with  time  period  i+1  ),  the  states  of  the  next 
stage  (associated  with  time  period  i  )  can  be  specified.  For  the 
first  stage  (associated  with  the  last  time  period),  since  there  is 
no  next  time  period,  the  states  are  not  constrained  b£  the  fourth 
type  of  constraint.  Therefore,  the  states  of  the  first  stage  of  the 
problem  are  given  as  a  result  of  finding  the  region  where  l  is 

the  last  time  period. 

For  this  problem  there  is  a  set  of  states  associated  with  the 
last  time  period.  Therefore,  it  is  necessary  to  determine  a  set  of 
optimal  policies,  each  associated  with  a  state  of  the  first  stage 
(last  time  period).  From  this  set  of  optimal  policies,  the  best  pol¬ 
icy  (lowest  cost)  is  the  true  optimum.  In  practice  this  is  accom¬ 
plished  by  solving  k  dynamic  programming  problems,  where  k  is 


i 

5 


j 


the  number  of  (c^.y^)  pairs  that  form  the  set  of  extreme  points 

of  the  region  that  satisfy  the  integer,  availability,  and  backward 
monotonicity  constraints  for  the  last  time  period. 

At  this  point,  all  of  the  techniques  needed  to  determine  the 
exact  solution  to  the  spares  provisioning  problem  have  been  dis¬ 
cussed.  These  can  be  summarized  as  follows. 

1.  Use  a  branch  and  bound  procedure  based  upon  the 

-  A  A  - 

approximation  of  X^  by  X^  (X^X^)  *  ^he 

A 

X^'s  are  independent  of  the  C  cj  >  Yj )  ,  J<i  , 
whereas  the  X^'s  are  not* 

2.  Enumerate  the  convex  hull  of  the  region  satisfy¬ 
ing  the  availability  constraint  and  integer  con¬ 
straint  for  each  time  period. 

3.  Use  corners  to  determine  the  D  regions.  This 
provides  a  method  of  enumerating  the  (c^,y^) 

pairs  associated  with  the  convex  hull  of  the  re¬ 
gion  satisfying  the  availability  constraint  for 
time  period  i  and  the  constraints  c^  >  ci_j  > 
y.  >  y.  .  . 

4.  By  dynamic  programming,  obtain  the  solution  to 
each  of  the  subproblems  associated  with  the  branch 
and  bound  procedure. 

The  next  chapter  will  discuss  the  implementation  of  these 
techniques  and  provide  computational  experience  associated  with 
implementation. 


CHAPTER  IV 


COMPUTATIONAL  EXPERIENCE 

This  chapter  discusses  the  applicability  and  utility  of 
the  approach  to  the  optimization  of  the  spares  provisioning  problem 
described  in  the  preceding  chapters.  All  of  the  techniques  dis¬ 
cussed  in  Chapter  III  have  been  incorporated  into  a  computer  program. 
A  listing  of  this  program  is  contained  in  Appendix  A.  After  a  dis¬ 
cussion  of  the  structure  of  the  program  and  the  function  of  its  ele¬ 
ments,  example  problems  are  presented.  The  nature  of  the  problems 
and  their  solution  are  then  discussed. 

The  computer  program  was  written  to  be  run  on  the  Hewlett 
Packard  3000  computer  of  the  School  of  Engineering  and  Applied  Sci¬ 
ence  of  The  George  Washington  University.  This  is  a  time-sharing 
system  that  allows  direct  input  of  data  at  the  time  the  program  is 
running.  This  latter  capability  is  utilized  in  performing  the 
branch  and  bound  procedure. 

The  program  consists  of  the  main  routine  and  three  subrou¬ 
tines.  The  main  routine  handles  the  data  input  and  the  output  of 
results,  as  well  as  the  logical  flow  of  the  solution.  Subroutine 
QUE  calculates  all  terms  associated  with  the  queueing  aspects  of  the 
problem.  This  includes  the  availability  constraint  and  the  mean 
number  repaired.  Subroutine  UNION  calculates  the  intersection  of 
dominated  regions  between  two  time  periods.  The  last  subroutine, 

DP,  generates  the  optimal  solution  of  a  branch  and  bound  subproblem 
through  dynamic  programming  techniques. 
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Inputs  to  the  program  for  each  time  period  are  the  failure 

rate,  A.  ;  the  mean  time  of  repair,  1/y .  ;  the  costs  C1  .  , 
i  X  X  j  1 

C_  ,  C  .  ,  and  C,  .  ;  and  the  availability  requirement  and  dis- 

count  factor.  As  discussed  in  Chapter  II,  only  the  costs  ^  , 

the  cost  associated  with  each  repair  channel,  and  ^  ,  the  cost 

associated  with  each  spare,  are  used  in  the  evaluation  of  the  objec¬ 
tive  function.  However,  for  each  solution,  the  actual  cost  based 

upon  C  .  ,  C  ,  C  ,  and  C,  .  is  also  calculated  and  printed 
i,i  i ,  i  -3,1  4,1 

out. 

A 

From  the  input,  the  A^'s  are  calculated.  These  are  lower 

bounds  on  the  A  's  ,  i.e.,  A,,  <  A .  ,  associated  with  the  extended 

i  1=1 

problem  of  the  branch  and  bound  procedure.  A  solution  is  obtained 

A 

using  the  A/s  and  then,  based  upon  the  optimal  policy,  the  cor¬ 
responding  A/s  are  calculated.  If  the  solution  satisfies  all  of 
the  constraints  when  they  are  evaluated  using  the  exact  A^s  >  then 

the  solution  is  feasible  and  the  problem  is  solved.  If  all  of  the 
availability  constraints  are  not  satisfied  using  the  exact  A^'s  » 

then  the  solution  is  infeasible  and  branching  must  therefore  occur. 

The  initial  experience  in  the  development  of  this  program 
indicated  that,  for  the  data  associated  with  the  gas  turbine  spares 
provisioning  problem,  the  solution  to  the  first  subproblem  of  the 
branch  and  bound  procedure  was  usually  feasible.  Therefore,  instead 
of  devising  a  complex  computer  program  that  would  generate  the  com¬ 
plete  branch  and  bound  tree,  it  was  decided  that  the  branching  pro¬ 
cedure  could  best  be  handled  external  to  the  computer,  using  the 
computer's  interactive  capability  to  input  the  parameters  associated 
with  each  new  branch.  A  later  modification  to  the  program  allows  the 
computation  of  the  solution  of  a  subproblem  and  the  set  of  branches 
that  directly  emanate  from  that  subproblera. 
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Because  the  complete  branch  and  bound  structure  is  not  com¬ 
puted  within  the  program,  the  Best  Lower  Bound  (BLB)  is  not  auto¬ 
matically  determined.  The  program  prints  out  the  optimal  value  of 
the  objective  function  and  the  value  of  the  total  cost,  as  well  as 
the  optimal  policy  for  each  of  the  branch  and  bound  problems.  It 
is  noted  if  the  solution  is  feasible.  The  (c^.y^)  Palrs  and  the 

associated  are  output.  The  program  does  keep  track  of  the 

Best  Upper  Bound  (BUB). 

Based  upon  the  output,  it  is  up  to  the  operator  to  deter¬ 
mine  the  BLB.  The  operator  must  then  input  the  next  problem  to  be 
solved.  This  is  accomplished  by  specifying  the  branching  time  per¬ 
iod  and  the  policy* for  the  previous  time  periods.  The  policy  for 
the  first  time  period  doe^  not  have  to  be  specified  because  the 
program  automatically  determines  it. 

Based  upon  this  new  input  problem,  the  program  evaluates 
a  solution,  determines  if  it  is  feasible,  branches  on  the  (c^.y^) 

pairs,  and  solves  the  new  problems  formed  by  the  branching.  If  the 
solution  to  a  problem  is  feasible,  no  further  branching  will  occur 
from  this  problem.  If  the  solution  to  a  feasible  problem  is  less 
than  the  BUB,  the  BUB  is  set  equal  to  the  value  of  the  objective 
function  for  this  new  feasible  solution.  If  the  value  ol  the  objec¬ 
tive  function  for  a  problem,  whether  it  is  feasible  oi  not,  is 
greater  than  the  BUB,  no  further  branching  need  occur  from  this 
problem.  This  is  because  the  theory  of  branch  and  bound  states  that 
neither  this  problem  nor  any  of  its  branches  can  yield  the  solution. 

The  procedure  is  best  illustrated  through  an  example.  Three 
test  problems  have  been  created.  Each  problem  has  five  time  periods 
with  identical  M^'s  ,  A^'s  ,  and  1/y^'s  ;  they  differ  only  in 

.  -  ■■•’its.  Table  IV-1  displays  the  values  used  for  each  , 

A  The  real  costs  and  the  discounted  costs  are  given 

•  '  the  computer  program,  the  A^’s  were  chosen  to 

r  ime  .is  opposed  to  most  actual  situations,  where 
■>pr  Ives  with  time. 

t 


i 

t 
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*  TABLE  IV-1 


INPUT 

PARAMETERS 
PROBLEMS  A, 

FOR  EXAMPLE 
B,  C 

Mi 

Ai 

Vv± 

10 

.00050000 

50.0000 

20 

.00060000 

50.0000 

30 

.00070000 

50.0000 

40 

.00070000 

50.0000 

50 

.00070000 

50.0000 

for  each  of  the  three  problems  in  Table  IV-2.  It  should  be  noted 


that  for  Problem  A,  C.  ,  =  CL  . 

1 ,  i  i ,  1 


J  ^  ;  for  Problem  B, 


.5C_  .  =  C  .  =  C  =  C  .  ;  and  for  Problem  C, 

^,1  J ,  1  111  4,1 


'5Cl,i  '  C2.1 


C3,i  C4,i 


The  computer  model  was  used  to  obtain  the  exact  so¬ 
lution  that  minimized  the  objective  function 


ci-l) 


+  C 


2,i 


(yi"yi-l 


This  is  labelled  OBJ  in  the  subsequent  figures.  The  total  cost, 

including  the  terms  C_  ,  and  C,  ,  is  labelled  COST  in  these 

J  ,  X  4  >  X 

figures . 


Figure  IV-1  shows  the  branch  and  bound  structure  associated 
with  Problem  A.  At  the  top  of  this  figure  is  the  extended  problem. 

A  A  - 

This  problem  is  solved  using  X^  ,  where  X^  <  X^  .  The  solution  is 
tested  for  feasibility  using  the  X^  calculated  from  the  specified 

A 

policy.  It  is  seen  that  this  solution  is  not  feasible.  Since  X^  «* 
A  -  A  ,  it  is  not  necessary  to  branch  at  the  first  time  period. 


For  this  problem  there  are  two  pairs,  (1,4)  and  (2,3).  These  pairs 
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TABLE  IV- 2 

COSTS  FOR  EXAMPLE  PROBLEMS  A,  B,  C 


Problem  A 

Costs 

c. 

r 

c„ 

c. 

1 

2 

3 

4 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

Discount 

Factor=.l 

Avail ability= . 9 

10.00 

10.00 

10.00 

10.00 

9.09 

9.09 

9.09 

9.09 

8.26 

8.26 

8.26 

8.26 

7.51 

7.51 

7.51 

7.51 

6.83 

6.83 

6.83 

6.83 

Problem  B 

Costs 

C, 

Cn 

C 

c. 

1 

2 

3 

4 

10.00 

20.00 

10.00 

10.00 

10.00 

20.00 

10.00 

10.00 

10.00 

20.00 

10.00 

10.00 

10.00 

20.00 

10.00 

10.00 

10.00 

20.00 

10.00 

10.00 

Discount 

Factor®. 1 

Aval lability®. 9 

10.00 

20.00 

10.00 

10.00 

9.09 

18.18 

9.09 

9.09 

8.26 

16.53 

8.26 

8.26 

7.51 

15.03 

7.51 

7.51 

6.83 

13.66 

6.83 

6.83 

Problem  C 

Costs 

c. 

c„ 

c„ 

c. 

1 

2 

3 

4 

20.00 

10.00 

10.00 

10.00 

20.00 

10.00 

10.00 

10.00 

20.00 

10.00 

10.00 

10.00 

20.00 

10.00 

10.00 

10.00 

20.00 

10.00 

10.00 

10.00 
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TABLE  IV- 2 — continued 


Discount 

Factor*. 1 

Availability*. 9 

20.00 

10.00 

10.00 

10.00 

18.18 

9.09 

9.09 

9.09 

16.53 

8.26 

8.26 

8.26 

15.03 

7.51 

7.51 

7.51 

13.66 

6.83 

6.83 

6.83 

are  associated  with  the  corners  defining  the  dominated  region  of  the 
second  time  period. 

Two  new  problems  are  created.  In  Problem  2,  the  solution 
must  contain  the  policy  =  1  ,  y^  =  4  •  The  exact  A^  is  calcu¬ 
lated  based  upon  this  policy.  The  subsequent  A^'s  (A^  and  A,.) 

A  A  — —  -A  - 

are  approximated  by  A^  such  that  A^  £  A^  and  A,.  <  A,.  .  For 

Problem  3  the  same  procedure  is  used  except  c?  =  2  and  y  ^  =  3  . 

-  A 

The  exact  A„  is  calculated  on  this  policy  and  A.'s  are  used  for 
3  x 

time  periods  4  and  5.  Problems  2  and  3  are  then  solved  and  tested 

for  feasibility. 

We  see  that  neither  Problem  2  nor  3  is  feasible.  Our  BLB  is 
70.79  and  there  is  no  BUB.  Problems  1,  2,  and  3  are  automatically 
solved  in  one  pass  of  the  computer  run. 

Since  the  BLB's  are  the  same  for  Problems  2  and  3,  but  the 
lower  cost  is  associated  with  Problem  2,  we  branch  on  this  problem. 
The  computer  program  allows  us  to  input  the  specified  policy  for 
time  periods  2  and  3.  It  is  not  necessary  again  to  input  the  values 
for  the  costs,  M^'s  ,  A^'s  ,  or  l/y^'s  . 

From  Problem  2,  in  time  period  2  there  are  three  possible 
policies,  (1,17),  (2,4),  and  (3,4).  The  computer  next  solves  Prob¬ 
lem  4,  which  was  created  by  inputing  C£ *  1  ,  y ^  »  c3  ■  1  •  and 


A  OBJ  70.79 
2  COST  375.73 


o  m 
n*  r-  b 
co  2 


>-i  cm  m  co  «<r  <r 


o  h  cm  co  co  <r 


O  m 
co  2 

>-i  cn  ro  <j-  <j-  ^3- 


U  H  CN  M  fO  vj 


o  m 

r-*  U> 
co  2 

>-<  cni  co  co  <r 


u  h  cm  cn  co  <r 


n  ro  vj  vj 


u  h  m  ro  CO  <f 


i-(  VO 

lu 
CO  2 

>♦  cm  <r  <r  <r  <r 


o  m  ^ 

Pn 

<n  2 

>-»  cvj  «**  ^  <r 


o  «h  h  co  co  <r 


o  m 
r>>  n* 
ro  2 


\0  r|  H  M  O 


. — ©: 


>h  cn  <r  <r  vj  <r 


u  h  hn  Osr 


o  in  y 

^  r-  r  P 

<n  fin  g 

n  sr  ^  m  ^ 
cn 

U  H  H  CM  co  CO 


00  CSJ 

co 

,  CO  fJH 

><  N  vO  vO 


CO 

vO  sO 
rH 

N  N  N  N 


U  H  H  fM  (N  fO 


U  rt  H  H  N  (SI 


Extended  branch  and  bound  tree  for  Problem 
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y^=17  .  Based  upon  this  policy  A ^  is  calculated  exactly  and  A^ 

is  determined.  The  computer  solves  this  problem.  Since  it  is  fea¬ 
sible,  no  further  branching  is  necessary. 

Problem  5  is  then  input  to  the  computer.  For  this  problem, 
c2 = 1  ,  y2 = 4  ,  c3 = 2  ,  and  y^ = 4  are  specified.  The  exact  A^ 

A 

is  calculated  and  A,,  determined  for  this  policy.  A  nonfeasible 

solution  is  found  and  two  more  branches  are  automatically  generated. 
We  see  that  the  solutions  to  Problems  6  and  7  are  both  feasible. 

For  Problem  7  the  value  of  the  objective  function  is  equal  to  the 
BLB.  Since  this  is  a  feasible  solution  we  know  that  it  is  an  exact 
solution  to  the  original  problem. 

At  this  point  the  branch  and  bound  procedure  could  be  termi¬ 
nated.  There  may  be  other  solutions  that  are  as  good  as  Problem  7, 
but  there  will  be  none  better.  For  the  sake  of  completeness,  the 
procedure  is  continued.  Problems  8,  9,  and  10  are  input  sequential¬ 
ly  to  the  computer.  Since  none  of  them  gives  feasible  solutions 
that  are  better  than  that  of  Problem  7  (which  determined  the  BUB) , 
no  further  branching  occurs. 

We  therefore  have  an  exact  optimal  solution  to  Problem  A 
given  by  Problem  7.  The  minimum  value  of  the  objective  function  is 
70.79  and  the  cost  associated  with  this  solution  is  375.51. 

Since  this  procedure  attempts  to  minimize  the  objective  func¬ 
tion  which  is  associated  with  only  the  direct  cost  of  the  repair 

channels  and  spares  (terms  associated  with  the  c.  .  and  c~  .),  it 

l,i  L , 1 

is  not  necessary  that  the  total  cost  (including  c_  R  )  associated 

J  >  i  i 

with  the  solution  be  the  minimum  possible  total  cost.  For  this 
problem  the  solution  provides  not  only  the  minimum  of  the  objective 
function,  but  the  minimum  total  cost  for  all  of  the  subproblems 
calculated. 

ProbLem  B  is  similar  to  Problem  A  except  that  the  cost  of 
the  spares  lias  been  doubled.  A  branching  procedure  similar  to  that 
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of  the  previous  problem  is  again  followed.  As  seen  in  Figure  IV-2, 
at  the  end  of  the  first  two  stages  we  have  no  feasible  solutions  and 
the  BLB  is  107.29.  Since  the  BLB  is  associated  with  Problem  3  we 
branch  on  iu.  When  Problem  4  is  input,  the  program  automatically 
generates  the  solutions  to  Problems  4,  5,  and  6.  We  next  input 
Problem  7,  generating  solutions  to  Problems  7  and  8.  This  exhausts 
the  branches  from  Problem  3. 

We  next  branch  on  Problem  2  by  inputting  Problem  9.  This 
has  a  feasible  solution  so  no  branching  comes  from  it.  We  then 
input  Problem  10  which  branches  into  Problems  11  and  12.  Finally, 
we  input  Problem  13  with  its  branch.  Problem  14. 

After  all  14  problems  have  been  solved,  each  of  the  solu¬ 
tions  to  the  unbranched  problems  is  feasible.  Therefore,  the  BLB 
equals  the  BUB  and  the  solution  is  associated  with  the  problem  that 
gives  the  BU(L)B.  Problem  8  gives  the  best  solution  and  again, 
serendipitously,  it  also  has  the  lowest  total  cost. 

From  Figure  IV-3  we  see  that  Problem  C  requires  the  solution 
of  only  one  problem.  For  this  case  the  solution  to  the  extended 
problem  is  feasible  and  therefore  an  exact  optimal  to  the  given 
problem. 

It  is  of  interest  to  compare  the  solutions  to  the  three 
problems.  Since  all  three  problems  use  the  same  values  for  the 
Mi's  ,  A.'s  ,  and  1/p^'s  ,  the  feasible  regions  are  identical. 

Therefore  the  enumeration  of  potential  solution  policies  is  the  same 
for  each  problem.  The  only  reason  the  solutions  will  differ  is  be¬ 
cause  of  the  differences  in  the  objective  function. 

Comparing  the  solutions  for  Problem  A  and  Problem  B  (Figures 
IV-1  and  IV-2,  respectively),  we  see  that  the  policies  for  the  first 
three  phases  of  the  branch  and  bound  are  identical.  The  costs  and 
values  of  the  objective  function  are  different.  Also,  if  a  policy 


OBJ  107.39 
COST  412.33 


Extended  branch  and  bound  tree  of  Problem 


C  Y  OBJECTIVE  96.57 

1  2  COST  403.74 

1  4  F 

2  4 

3  4 

3  5 

Fig.  IV-3.  Branch  and  bound  tree  for  Problem  C. 


is  not  feasible,  it  is  not  feasible  for  all  three  problems  (A,  B, 
and  C) .  Problems  A  and  B  differ  only  in  the  last  phase  of  the 
branch  and  bound  structure.  In  Problem  A,  because  of  the  order  in 
which  the  branching  was  taken,  only  ten  problems  were  solved.  In 
fact,  since  a  feasible  solution  was  found  at  Problem  7  whose  value 
of  the  objective  function  was  equal  to  the  BLB,  only  seven  problems 
needed  to  be  solved.  For  Problem  B,  since  no  feasible  solution  had 
the  value  of  the  objective  function  equal  to  the  BLB  of  the  previous 
phase,  all  branches  had  to  be  evaluated. 

By  looking  at  the  three  problems  (A,  B,  and  C) ,  we  see  that 

it  is  not  obvious  how  efficient  the  branch  and  bound  procedures  may 

be.  For  Problem  C  it  was  necessary  only  to  solve  one  problem.  For 

Problem  B  all  of  the  branches  had  to  be  calculated.  Problem  A  could 

be  solved  with  only  about  half  of  the  branches  evaluated.  While  it 

is  impossible  to  predict  exactly,  it  is  suspected  that  the  closer 

the  bounds  on  A .  ,  i.e.,  A.  ,{M,  -./M.lfA.  -  A.  ,l«0  ,  for 

l  i-l>  i-1  l*'  i-1  l-l ' 

>.  ,  or  i  ~  ^  j)  ft*  0  for  [Equation 

(111-1)1,  the  earlier  the  branch  and  bound  procedure  will  obtain  the 
solution. 

The  primary  motivation  for  this  dissertation  is  based  upon 
the  work  given  in  Reference  [1].  In  that  paper  a  specific  gas  tur¬ 
bine  problem  was  solved.  The  input  parameters  associated  with  this 
problem  are  given  in  Table  IV-3.  This  is  an  eleven  time-period 
problem;  all  costs  are  given  in  $1,000  units. 
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TABLE  IV- 3 

INPUT  TO  GAS  TURBINE  PROBLEM 


M. 

i 

X. 

l 

1/Pi 

10 

.00147186 

65.0000 

28 

.00152455 

62.5000 

50 

.00136767 

60.0000 

82 

.00099101 

57.5000 

121 

.00082569 

55.0000 

158 

.00071831 

55.0000 

182 

.00063699 

55.0000 

208 

.00060941 

55.0000 

229 

.00061725 

55.0000 

251 

.00062297 

55.0000 

256 

.00062015 

55.0000 

Costs 

C, 

C„ 

C„ 

c. 

1 

2 

3 

4 

132.00 

822.00 

49.00 

1975.00 

132.00 

945.00 

49.00 

2760.00 

132.00 

1087.00 

37.00 

3840.00 

132.00 

1174.00 

40.00 

3950.00 

132.00 

1268.00 

42.00 

2800.00 

132.00 

1369.00 

44.00 

1100.00 

132.00 

1369.00 

44.00 

350.00 

132.00 

1369.00 

44.00 

350.00 

132.00 

1369.00 

44.00 

350.00 

132.00 

1369.00 

44.00 

350.00 

132.00 

1369.00 

44.00 

350.00 

Discounted  Costs:  Factor®. 1 

C, 

C„ 

C„ 

c. 

1 

2 

3 

4 

132.00 

822.00 

49.00 

1975.00 

120.00 

859.09 

44.55 

2509.09 

109.09 

898.35 

31.24 

3173.55 

99.17 

882.04 

30.05 

2967.69 

90.16 

866.06 

28.69 

1912.44 

81.96 

850.04 

27.32 

683.01 

74.51 

772.77 

24.84 

197.57 

67.74 

702.51 

22.58 

179.61 

61.58 

638.65 

20.53 

163.28 

55.98 

580.59 

18.66 

148.43 

50.89 

527.81 

16.96 

134.94 

Availability® ,9 


»■:  'vr’«S**5!h> '*•- JUS.--** 
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Using  the  exact  solution  computer  program,  the  solution  at 
the  top  of  Table  IV-4  was  obtained.  This  is  compared,  in  this 
table,  to  the  solution  generated  using  the  heuristic  approach  of 
Reference  [1],  We  see  that  the  exact  solution  is  slightly  more 
than  $300,000  less  costly.  The  major  difference  between  the  optimal 
policies  is  that  the  exact  solution  buys  spares  in  an  earlier  time 
period,  when  they  are  cheaper,  rather  than  waiting  until  they  are 
needed  in  a  later  time  period.  In  addition,  in  the  heuristic  solu¬ 
tion,  repair  channels  are  purchased,  given  up,  and  then  repurchased 
(time  periods  8,  9,  10). 

The  procedure  itemized  in  this  paper  is  not  only  exact  but 
also  efficient.  For  the  gas  turbine  problem  only  one  problem  had 
to  be  solved.  There  was  no  branching  necessary.  A  total  of  only 
66  seconds  of  Hewlett  Packard  3000  CPU  time  and  6  minutes  of  termi¬ 
nal  time,  including  input  of  the  data  and  output  of  the  results, 
was  required  to  solve  the  problem.  This  procedure  has  been  tested 
on  only  a  limited  number  of  problems,  but  it  is  expected  that  it 
will  be  significantly  more  efficient  than  the  heuristic  procedure 
where  more  interaction  and  a  greater  number  of  iterations  are  re¬ 
quired.  How  this  computational  efficiency  can  be  extended  to  other 
problems  will  be  discussed  in  the  next  chapter. 
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TABLE  IV- 4 

SOLUTION  TO  GAS  TURBINE  PROBLEM 


COST  =*  38827.16 


Icy  R 


1 

2 

8 

5.371 

2 

4 

8 

15.337 

3 

8 

8 

26.426 

4 

10 

10 

37.197 

5 

10 

12 

45.492 

6 

12 

13 

51.798 

7 

12 

14 

53.967 

8 

12 

14 

56.266 

9 

13 

14 

58.288 

10 

15 

14 

61.583 

11 

15 

14 

61.600 

Exact  Solution 


c 

y 

R 

Present 

Worth 

3 

3 

5.3 

5097.61 

5 

6 

15.3 

11105.17 

8 

8 

26.4 

17226.10 

10 

10 

37.2 

23273.15 

12 

11 

45.5 

27537.26 

15 

12 

51.8 

30731.45 

13 

13 

53.9 

33041.68 

14 

13 

56.3 

34559.60 

13 

14 

58.3 

36557.30 

15 

14 

61.6 

37966.66 

15 

14 

61.6 

39146.38 

Heuristic  Solution 


CHAPTER  V 


SUMMARY,  PROBLEM  EXTENSIONS ,  AND  CONCLUSIONS 


Summarv 


A  methodology  for  attaining  an  exact  solution  to  the  problem 
of  minimizing  the  cost  of  a  multi-year  spares  provisioning  problem 
has  been  developed  in  this  dissertation.  The  problem  solved  is  a 
nonlinear  integer  programming  problem,  and  was  previously  solved 
using  a  heuristic  algorithm.  There  was  no  assurance  that  the  heu¬ 
ristic  solution  was  optimal  or  even  near  optimal.  However,  compari¬ 
son  of  results  of  this  effort  and  the  heuristic  indicate  that  the 
heuristic,  in  the  cases  tested,  gives  results  very  close  to  optimal. 

There  are  a  number  of  difficulties  associated  with  this 
problem,  beyond  those  normally  associated  with  integer  programming 
problems.  The  forms  of  the  equations  of  state  are  dependent  not 
only  on  the  state  variables,  but  also  on  the  relative  values  of 
these  variables  as  well  as  other  parameters  of  the  problem.  In 
addition,  a  term,  X  ,  which  is  normally  a  parameter  of  the  problem, 
is  a  function  of  the  state  variables.  To  resolve  these  problems, 
an  amalgamation  of  techniques  has  been  used  to  obtain  the  solution. 
Branch  and  bound  is  used  to  reparameterize  X  .  Enumeration  and  a 
construct  of  "corners"  are  used  to  specify  the  feasible  region  of 
the  subproblems.  Finally,  dynamic  programming  is  used  to  obtain  a 
solution.  This  dissertation  demonstrates  that  this  approach  gives 
an  exact  optimal  solution  and  that  this  solution  is  generated  in  an 
efficient  manner. 
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This  effort  was  motivated  by  the  need  to  obtain  an  exact 
solution  to  a  specific  problem,  the  minimum  cost,  multi-year 
policy  for  purchasing  spares  and  repair  channels  for  a  gas  turbine 
maintenance  system.  There  is  no  reason  why  the  applicability  of 
the  solution  technique  is  limited  to  the  particular  gas  turbine 
problem.  The  underlying  stochastic  model,  the  machine  repair 
problem  with  spares,  is  the  foundation  for  analysis  of  many  sys¬ 
tems  problems.  This  chapter  discusses  the  types  of  problem  exten¬ 
sions  that  can  be  handled  with  minor  modifications  to  the  solution 
framework.  Also  discussed  is  the  potential  applicability  of  some 
of  the  concepts  developed  in  obtaining  the  solution  techniques. 

V. 2 _ Problem  Extension 


There  are  three  basic  factors  that  determine  the  nature  of 
the  solution  to  this  problem.  They  are: 

1.  a  linear  objective  function; 

2.  a  stochastic  availability  constraint  based  upon 
the  equations  of  state  associated  with  the  clas¬ 
sical  machine  repair  problem  with  spares; 

3.  a  key  parameter  in  determining  the  equations  of 
state,  A^  ,  the  mean  failure  rate,  which  is  an 
implicit  function  of  the  decision  variables. 

Each  of  these  factors  had  significant  impact  on  the  choice  of 
techniques  required  for  an  exact  solution.  The  linearity  of  the  ob¬ 
jective  function  forces  the  solution  to  lie  on  the  convex  hull  of 
the  feasible  region.  This  allows  an  enumeration  of  possible  solution 
sets  because  the  region  can  be  described  by  a  finite  number  of 
points.  The  use  of  the  solution  to  the  machine  repair  problem  pro¬ 
vided  a  functional  dependence  on  the  state  variables  that  could  be 
enumerated.  Finally,  the  dependence  of  A^  on  the  state  variables 

requires  the  use  of  the  branch  and  bound  procedure  to  approximate 
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successively  the  problem  in  a  manner  that  assures  an  exact  solu¬ 
tion. 

Dynamic  programming  is  not  necessary  for  obtaining  a  solu¬ 
tion.  Once  the  convex  hull  of  the  feasible  region  has  been  enu¬ 
merated,  it  could  be  described  by  linear  inequalities,  allowing 
solution  through  linear  programming  techniques.  However,  dynamic 
programming  provides  computational  efficiency  and  latitude  in  modi¬ 
fication  of  the  types  of  problems  that  can  be  solved. 

The  first  possibility  is  modification  of  the  objective 
function.  Its  linearity  was  used  only  to  guarantee  that  the  prob¬ 
lem's  solution  is  on  the  convex  hull.  Since  any  objective  function 
that  is  isotonic  increasing  [a  function  is  isotonic  if,  for  vectors 
a  and  b  where  a^>b  ,  then  f(a)^f(b)]  in  the  state  variables 
will  also  guarantee  this,  such  extended  problems  can  be  solved 
through  the  procedure  described. 
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C  .  *  operating  cost  of  a  server,  year  i 
°,i 

C  *  inventory  carrying  cost,  year  i 

1,1 

C  .  •  cost  of  CIP,  year  i 
r  ,i 

(a)+  *  a  if  a>0  and  0  if  a^O  ; 

(a)  *  a  if  a<0  and  0  if  a_>0  . 

It  has  not  been  determined  if  this  objective  function  is  isotonic. 

It  is  expected  that  an  objective  function  containing  only  the  first 
five  terms  can  be  handled  using  the  technique  of  the  exact  solution. 
While  such  an  objective  function  is  not  necessarily  isotonic,  it 
should  be  possible  to  enumerate  all  potential  interior  solution 
points.  These  states  can  be  added  to  the  state  space  of  the  dynamic 
program.  Whether  the  entire  objective  function  can  be  handled  de¬ 
pends  on  the  nature  of  the  repair,  inventory  and  CIP  costs,  which 
requires  further  analysis.  If  such  an  objective  function  is  used, 
an  additional  modification  to  the  problem  would  have  to  be  made. 

The  constraints  c^  c^_^  and  y^  £  y^_^  are  no  longer  valid. 

The  generation  of  the  feasible  space  would  have  to  be  changed.  It 
is  expected,  however,  that  these  changes  could  be  made  within  the 
general  framework  of  the  solution  with  little  modification  of  the 
computer  program. 


Another  extension  of  the  original  gas  turbine  problem  is  a 
modification  of  the  definition  of  the  availability  constraint. 
Instead  of  requiring  a  minimum  spare  availability  at  the  time  of 
failure,  the  constraint  is  placed  on  the  overall  minimum  fleet 
availability  over  the  total  time  period.  This  is  equivalent  to 
replacing  the  constraints 


by 


y-1 


n**0 


n,  i 


a 


nm0 


Pn,i  - 


♦ 


where 


=  Pr{n  in  system | arrival  is  about  to  occur}  , 

p  =  Pr{n  in  system}  , 
n 

a  is  the  availability  required. 

This  new  constraint  is  computationally  simpler  to  work  with  and 
should  satisfy  the  important  property  that  if  the  availability  con¬ 
straint  is  satisfied  for  a  pair  (c,y)  then  it  is  satisfied  for  a 
pair  (c',y')  ,  where  c'^c  and  y'^y  . 

These  have  been  explicit  modifications  to  the  problem  that 
have  application  to  ongoing  work.  The  general  nature  of  the  problem 
with  spares  arises  in  industries  as  diverse  as  utilities  to  automo¬ 
bile  repair.  While  each  problem  has  its  own  features,  the  generality 
of  the  approach  should  allow  solution  of  a  number  of  different  types 
of  problems. 

Another  aspect  of  this  dissertation, that  may  have  use  beyond 
this  immediate  application,  is  the  concept  of  "corners."  In  this 
dissertation  very  simple  relationships  between  point  sets  were 
handled  in  a  computationally  efficient  manner.  It  may  be  possible 
to  handle  more  complex  relationships  with  the  same  efficiency.  This 
technique  simplifies  the  enumeration  of  states  associated  with  the 
integer  programs.  Further  exploration  of  this  concept  may  expedite 
the  solution  of  other  types  of  integer  programs. 

V.3  Conclusion 


In  an  attempt  to  obtain  an  exact  solution  to  a  nonlinear 
integer  programming  problem,  a  diverse  group  of  techniques  was 
brought  together  into  a  single  package  that  is  easy  to  use  and 
computationally  efficient.  It  is  hoped  that  these  techniques  can 
find  wider  application  in  the  general  area  of  problems  associated 
with  the  classical  problem  of  machine  repair  with  spares. 


APPENDIX  A 


LISTING  OF  COMPUTER  PROGRAM 


90 


1  C  THIS  PROGRAM  DETERMINES  THE  OPTIMAL  POLICY  FOR  A  SPARES 

2  C  PROVISIONING  PROBLEM.  A  BRANCH  AND  BOUND  APPROACH  IS 

3  C  USED.  THEREFORE  SUBPROBLEMS  MAY  HAVE  TO  BE  INPUT  SERUENT- 

4  C  IALLY .  THIS  APPROACH  ALSO  INCORPORATES  ENUMERATION 

5  C  TECHNIQUES  AND  DYNAMIC  PROGRAMING. 

6  c************************************************************* 

7  c  input: 

8  C  NT=NUMBEROF  TIME  PERIODS 

9  C  M( I) -NUMBER  OF  MACHINES  *  TIME  PERIOD  I 

10  C  X ( I ) -LAMBDA  (FAILURE  RATE)  FOR  TIME  PERIOD  I 

11  C  XI ( I) -TIME  BETWEEN  REPAIRS  (1/MU)  FOR  TIME  PERIOD  I 

12  C  SC(I»J)=COSTS  FOR  TIME  PERIOD  I . J- 1 » 2 » 3 * 4 

13  C  ASSOCIATED  WITH  Cl t C2 > C3 » C4 » OF  THE  OBJECTIVE 

14  C  FUNCTION  RESPECTIVELY 

15  C  DIS-DISCOUNT  FACTOR 

16  C  AV=AVA I LABILITY  REQUIREMENT 

17  C  IP+PRINT  SWITCH f -1  FULL  PRINT f  OTHERWISE  PRINT  SOLUTION 

18  C  OUTPUT  ONLY 

19  C  IBB-BRANCHXNG  SWITCH f -0  STOP i OTHERWISE  THE  NUMBER  OF 

20  C  THE  BRANCHING  TIME  PERIOD 

21  C  IX(Krlrl)-NUMBER  OF  REPAIR  CHANNELS  SPECIFIED  FOR  TIME 

22  C  PERIOD  K  FOR  BRANCHING  PROBLEM. 

23  C  1<K<IBB 

24  C  I X ( K » 1 » 2 ) -NUMBER  OF  SPARES  SPECIFIED  FOR  TIME  PERIOD 

25  C  K  FOR  BRANCH  PROBLEM. 

26  C  1 <N<IBB 

27  C  R ( K ) “MEAN  NUMBER  OF  MACHINES  REPAIRED  IN  TIME  PERIOD 

28  C  K  FOR  BRANCHING  PROBLEM. 

29  C  1<K<IBB 

30  C******************************************** 

31  C  OUTPUT 

32  C  SC ( I > J ) ~C 1 » C2 » C3 » C4  FOR  J+lr2r3r4r  RESPECTIVELY »  FOR  TIME 

33  C  PERIOD  I 

34  C  -DISCOUNTED  COST  C1»C2»C3»C4 

35  C  -DIFFERENCE  IN  COST  BETWEEN  TIME  PERIODS  I  AND 

36  C  I  f  1 »  FOR  Cl  AND  C2»  J=1 » 2* RESPECTIVELY 

37  C**SC(I»J>  IS  USED  FOR  THREE  DIFFERENT  COSTS  TO  REDUCE 

38  C#*MEMMQRY  REQUIREMENTS  OF  THE  COMPUTER  PROGRAM** 

39  C  I X < I »  J » K ) -CORNERS  OF  THE  DOMINATED  REGIONS (SET  OF 

40  C  MINIMUM  PAIRS (C » Y )  SATISFYING  THE  CONSTRAINTS) 

41  C  T~1  TO  NT t  NUMBER  OF  TIME  PERIODS 

42  C  J-NUMBEE  OF  PAIRS  IN  THE  SET 

43  C  K=1 5 NUMBER  OF  REPAIR  CHANNELS 

44  C  2  f  NUMBER  OF  SPARES 

45  C  ID( I ) -NUMBER  OF  PAIRS  TN  SET  FOR  TIME  PERIOD  If  J<ID(I) 

46  C  Cl -TOTAL  COST  OF  THE  SOLUTION  INCLUDING  Cl »C2»C3fC4 

47  C  CO-COST  OF  THE  OBJECTIVE  FUNCTI0W»C1  AND  C2  TERMS  ONLY 

48  C  I  -- 1 1 MF  PERIOD 

49  C  JS ( I » 1  )  -NUMBER  OF'  REPAIR  CHANNELS  AVAILABLE  IN  TIME  PERIO 

I 

jfSM  COPY  FUMlJiifi"  lv  uyi. 


r  {  "  " 

• 

• 

91  9 

51 

c 

OF  THE  SOLUTION  TO  THE  SURPR0BLEM 

52 

c 

JS ( I , 2 ) =NUMBER  OF  SPARES  AVAILABLE  IN  TIME  PERIOD  I  OF  THE 

53 

c 

SOLUTION  TO  THE  SUBPROBLEM 

54 

c 

R< I >=MEAN  NUMBER  OF  MACHINES  REPAIRED  IN  TIME  PERIOD  I 

55 

c 

OF  THE  SOLUTION  TO  THE  SUBPROBLEM 

56 

57 

58 

c 

Z( I )=LAMBDA  BAR  (AVERAGE  FAILURE  RATE)  FOR  THE  TIME  PERIOD 

c 

I  ASSOCIATED  WITH  THE  SOLUTION  TO  THE  SUBPROBLEM 

c 

IBD=BRANCHING  TIME  PERIOD 

59 

c 

ICY=NUMBER  OF  REPAIR  CHANNELS 

60 

c 

IYY=NUMBER  OF  SPARES 

1  61 

c 

RT=MEAN  NUMBER  OF  MACHINES  REPAIRED 

62 

A 

O  ^  di  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  di  ^  ^  ^  ^  lb  ^  ^  X  X  X  .L  vLr  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  ■ 

w  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ip  ^  ^  ^  ^  ^  Ip  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  Jp  ^  ^  ^  V 

63 

c 

1 

64 

c 

i 

65 

C 

THIS  IS  THE  MAIN  PROGRAM.  IT  ACCEPTS  INPUT,  WRITES  OUTPUT 

66 

C 

, AND  GOVERNS  THE  FLOW  OF  THE  SOLUTION.  IT  AUTOMATICALLY 

67 

C 

GENERATES  ONE  SET  OF  BRANCHES  IF  REQUIRED. 

100 

DIMENSION  M(15),X(15)»Xl(i5),R(15),Z(15),JS(15,2), 

110 

1ID< 15), SC (15, 4), IX (15, 25, 2) ,R3(25) 

120 

WRITE (6, 910) 

130 

ACCEPT  NT 

140 

WRITE (6, 911) 

150 

DO  91  1=1, NT 

160 

91  ACCEPT  M(I),X(I),X1(I)  ! 

170 

Z(1)=X(1)  j 

180 

WRITE( 6 » 902 )  (M( I ) , X ( I ) , XI ( I ) , 1  =  1 ,NT)  j 

190 

WRITE (6,912)  \ 

200 

DO  92  1=1, NT  ! 

210 

92  ACCEPT (SC(I , J) , J=1 ,4) 

1  220 

WRITE (6,917) ( (SC(I,J) ,J=1 ,4) ,1=1, NT) 

230 

WRITE(6,913)  ! 

240 

ACCEPT  DIS, AV 

250 

WRITE (6, 918)  DIS , AV 

260 

DD-1 

270 

DO  1  1=2, NT 

280 

DD=DD/(DIS+1.)  j 

290 

DO  1  J=l,4  j 

300 

1  SC(I,J)=SC(I,J)*DD  ■■ 

310 

WRITE(6,917 )  ((SC(I,K),K=1,4),I=1,NT) 

320 

R  ( 1  )  =0 

325 

C 

CALCULATE  COST  DIFFERENCES-USED  IN  MODIFIED  OBJECTIVE  i 

326 

C 

FUNCTION  i 

330 

DO  2  1=2, NT 

340 

J=I-1 

350 

SC(J,1)=SC(J,1)-SC(I,.1.) 

360 

2  SC ( J , 2 ) =SC ( J , 2 ) -SC (1,2)  1 

370 

WRITE (6, 914) 

380 

ACCEPT  IP 

»  390 

a 

IF(IP.EQ.l)  WRITE ( 6 , 906 )  ( ( SC ( T , K ) , K~ t , 4 ) , I - 1 » NT ) 

I*  * 

,  ,  , 

t 

L* 

92 


400 
410 
420 
425 
430 
440 
450 
4  60 
470 
480 
490 
500 
505 
510 
520 
530 
540 
550 
555 
540 
570 
580 


BUB= 100000000 
IBB=1 
3  LS=1 

C  BRANCH  AND  BOUND  LOOP 
DO  250  IN=:l  »2 
DO  200  N=1 r LS 
IR=0 

DO  .100  K-l  »NT 

IF(K  « LE 4 1 )  GO  TO  4 

J=K-1 

I=K 

IFfMCD.LT.MfJ))  I=J 
C  CALCULATE  LAMBDA  CAP 

Z(K)-( <M(I)-M<J> )*X(K)+R( J>*X(J)+(M( J)-R( J) >*Z< J) >/M(I) 
4  CONTINUE 

IF(K,LT4lBB)  GO  TO  100 
R ( K ) =0 

I F ( Z < K ) .GE.X(K) >  R<K>=365.*Z(K)*M<K> 

C  CALCULATE  RHO 

XM=Z(K>*M<K>*X1<K> 

L=0 
I  Y0=0 


590 

600 

610 

620 

630 

635 

640 

650 

660 

665 

670 

680 

690 

700 

705 


IC0=0 

ICY=XM-1 

IF ( ICY  4  LT  4  0 ) ICY~0 
5  ICY~ICY+1 
I YY-ICY 

C  CALCULATE  MINIMUM  <C,Y>  PAIRS 

10  CALL  GUE ( XM » M  <  K ) f  AU » I C Y » I YY » I R » X 1 ( K ) »  RT  > 
IF < IP 4 EG. 1 )  WRITE  <  A » 901 )  ICY.IYY 
IF ( ICY , GT .  I YY+M < K ) )  GO  TO  70 
C  DETERMINE  IF  C  AND  Y  ARE  BOTH  NEW 

IF( (IYY~IYO>#< ICY- I CO) .EQ.O)  GO  TO  60 
L~L+1 
ICO= ICY 
IYO-=I  YY 

C  IX  IS  A  CORNER  POINT 


710 

720 

730 

735 

740 

745 

744 

750 

740 

770 

780 

790 

800 


IX(K>L» 1 >-ICY 
IX(K»L»2)=IYY 
60  CONTINUE 
C  TEST  IF  Y>C 

IF  (IYY.GT.ICY)  GO  TO  5 

C  TEST  IF  ALL  CORNER  POINTS  HAUE  BEEN  GENERATED  FOR  THIS 
C  TIME  PERIOD 

IF C IYY . EG. 1 >  GO  TO  70 
I YY~IYY-i 
GO  TO  10 
70  ID(K)~L 

IF ( IP . EG . 1 )  WRITE <6» 905  >  (K » IX< K» I » 1 ) » IX<K » I » 2 ) » 1  =  1 »L > 
100  CONTINUE 


IHIS  PAdS  IS  BEST  QUALITY  PRACTICABLE 

ffcOM  .X.rY  FURBISH  TO  DDC 
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810 

81S 

820 

830 

K= 1 » NT ) 
835 
840 
850 
855 
860 
870 
880 
890 
900 
910 
920 
925 
930 
940 
950 
960 
970 
975 
980 
990 
1000 
1010 
1020 
1030 
1040 
1045 
1050 
1060 
1070 
1075 
1080 
1085 
1090 
1100 
1105 
1110 
1120 
1130 
1140 
1145 


WRITE(6»903>  (Z( I ) > 1-1 »NT) 

C  MERGE  DOMINATED  REGIONS 
CALL  UNION (NT  > IX » ID ) 

IF(IP«EQ.1)WRITE(6»  905 ) ( <K»  < IX <K , I . J) . J=1 »2 ) » 1=1 » ID(K ) ) f 

C  DO  DYNAMIC  PROGRAMMING 

CALL  DP (NT» IX » IB* SC»  JS  » CS » IP) 

CT=CS 

C  TEST  IF  SOLUTION  IS  FEASIBLE  FOR  ORIGINAL  PROBLEM 
X0=0. 

RO-O . 

M0=0 

Z0=0. 

IRS-0 

IR=1 

DO  150  1=1 *NT 
C  CALCULATE  LAMBDA  BAR 

Z ( I ) * < ( M  < I ) -MO ) *X ( I ) +R0*X0+  < MO-RO ) *Z0 ) /M ( I ) 
XM=Z(I)*M(I)*X1(I> 

ICY=JS( 1*1) 

IYY=JS< 1*2) 

IR=1 

C  TEST  IF  AVAILABILITY  CONSTRAINT  IS  SATISFIED 
CALL  QUE(XM»M< I ) *  AV» ICY »IYY*IR»X1(I)»RT) 

R(I  )=365.JKRT 
XO=X(I) 

RO=R ( I ) 

MO=M ( I ) 

ZO=Z(I> 

IF  ( IR . LT . 0 )  IRS52! 

C  CALCULATE  TOTAL  COST 

150  CT---CT+R  ( I )  *SC  (1*3)  +SC  ( 1 * 4  ) 

WRITE <6* 908)  CT»C5» ( I  * JS( I  *  1 ) » JS< I»2)»R(I)*I-1*NT) 

WRITE ( 6 » 903 )  (Z( I ) * 1=1 *NT) 

C  TEST  IF  FEASIBLE.  IF  INFEASIBLE  BRANCHING  IS  REQUIRED 
IF < IRS.EQ.O)  WRITE (6.916) 

C  TEST  IF  THE  SOLUTION  IS  A  BEST  UPPER  BOUND 
IF(CT.GT.BUB)  GO  TO  170 
IF(IRS.NE.O)  GO  TO  155 
C  HAS  BRANCHING  ALREADY  OCCURED 
IF ( IBB ♦ EG • 1 >  STOP 
BUB»CT 
GO  TO  170 
155  CONTINUE 
C  BRANCH 


_  T  r-  rtE3T  auftlitl  metiers 
oot  *  *  *v 
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1  150 
1155 
1160 
1170 
1180 
1185 
1190 

1195 

1196 
1200 
1210 
1220 
1230 
1235 
1240 
1250 
1260 
1270 
1280 
1290 
1295 
1300 
1310 
1320 
1330 
1335 
1340 
1350 
1360 
1370 
1380 
1385 
1390 
1400 
1410 
1420 
1430 
1440 
1445 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 


IF( IBB.EQ. 1 )  TBB-2 

C  DETERMINE  TIME  PERIOD  FOR  BRANCHING 
IBD-IBB+l-IN 
WRITE < 6 »  909 )  IBD 
L=ID<IBB> 

C  CALCULATE  RHO 

XM~Z ( IBB  >#M  < IBB) #X1 ( IBB  > 

C  DETERMINE  LAMBDA  BAR  FOR  THE  TIME  PERIOD  PRIOR  TO  THE 

C  BRANCHING  TIME  PERIOD 
DO  160  I  ~  1  *  L 
IC.Y~IX(  IBD » I » 1  > 

I YY~IX ( IBB » I » 2 ) 

IR~1 

C  CALCULATE  LAMBDA  BAR 

CALL  QUE <  XM  t M < IBB >  » AV » ICY » I YY r IR  r  XI ( IBB  > > RT ) 

RT=RT*365 

IF < IN . EQ ♦ 1 )  RS(I)=RT 

160  WRITE <6> 915)  ICY» IYY»RT 

170  LN=IN+N-1 

IF(LN.GT.LS)  GO  TO  200 

C  DETERMINE  BRANCHES 

I X ( IBD > 1 * 1 )~IX( I BD > LW » 1 ) 

IX  < IBD » 1 » 2  >  *IX  < IBD » LN » 2 ) 

R<IBD>=RS<LN> 

200  CONTINUE 

C  TEST  IF  THIS  BRANCH  IS  GREATER  THAN  THE  BUB 
IFCCT.GE.BUB)  GO  TO  300 
ID< IBB ) ~ 1 
::  BB-iBB+i 

250  LS-L 

300  WRITE ( 6  r 904 ) 

C  INPUT  NEW  BRANCH 
ACCEPT  IBB 
IF ( IBB ♦ LE . 0 )  STOP 
IF <1  BB.L.T *3)  GO  TO  4 
J-IBB-t 
DO  35  K'-2 1 J 
WRITE  <  6  v  907 )  K 

C  INPUT  POLICY  ASSOCIATED  WITH  NEW  BRANCH 
ACCEPT  I X  <  K  V 1 y 1 ) fix ( K rl r  2) r R ( K ) 

35  I D  ( K  >  -- 1 
GO  TO  3 

?  0 1  F  0  R  M  A  T  ( 2 1 5  r  F 1 0 .  .3 ) 

902  FORMAT ('  M  LAMBDA  i/rUJ  VI 1  <  I6»F1 1 .8  »F8. 4/ >  } 

903  FORMAT  ( // /  LBAR  ■ '  » 6F:i  0  .S/6X ,  6F:I  0 .8//) 

904  FORMAT (//'  INPUT  BRANCHING  TIME  PERIOD ’ > 

905  FORMAT  (  '  K='  »  .T5r  '  C«  '  » 15 »  '  Y~  '  r  15  > 

906  FORMAT (  4F10.2) 

907  FORMAT <'  INPUT  CrY  AND  RfiAR  FOR  TIME  PERIOD  ’ » T ? ) 

900  FORMAT ( /3X r ' COST- ' »F10.2»3X*  'OBJECT TOE™ ' *F10.:V 

THIS  PASS  IS  BEST  QUALITY  PRACTICED* 
nouM  COII  10  DDC  ...  — ' 


ns 


-95 


1560 

1570 

1580 

1590 

1600 

1610 

EMENT'/) 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1681 

1682 

1683 

1684 

1685 

1686 

1687 

1688 

1689 

1690 
1700 
1705 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1785 
1790 
1800 
1810 
1815 
1820 
1830 
1840 
1850 
1855 
1860 
1870 
1880 
1890 
1900 
1910 


12X»'I'»4X>'C'»4X»'Y'»7X»'R'/<3I5»F10»3)) 

909  FORMAT ( '  BRANCHING  ON  TIME  PERIOD  '»I3> 

910  FORMAT <5Xt 'ENTER  NUMBER  OF  TIME  PERIODS'/) 

911  F0RMAT(5X» 'ENTER  M,X>1/MU'/) 

912  F0RMAT<5X> 'ENTER  COSTS'/) 

913  FORMAT <5X* 'ENTER  DISCOUNT  FACTOR  AND  AVAILABILITY  REQUIR 

914  FORMAT (5X » ' IF  DIAGNOSTIC  PRINT  DESIRED  ENTER  1') 

915  FORMAT ( '  C='»I4*'  Y=',I4,'  RBAR=  ' , F12 ♦ 8 ) 

916  FORMAT < '  FEASIBLE  SOLUTION'/) 

917  FORMAT <4X> 'Cl ' ,8X> 'C2' »8X> 'C3' »8X* 'C4'/15(4F10.2) > 

918  FORMAT ( '  DISCOUNT  FACTOR= ' » F6 ♦ 4  f '  AVAILABILITY- ' »F6 ♦ 4 ) 
END 

SUBROUTINE  UNION < NT > IX » ID) 

THIS  SUBROUTINE  CALCULATES  THE  CORNERS  GENERATED  BY 
TAKING  THE  INTERSECTION  OF  TWO  DOMINATED  REGIONS 
NT -THE  NUMBER  OF  TIME  PERIODS 
IX < I » J»K)-THE  < C » Y >  PAIR  DEFINING  THE  CORNERS 
I=THE  TIME  PERIOD 
J=THE  CORNER  INDEX 
K=1  FOR  C » 2  FOR  Y 

ID< I )=THE  NUMBER  OF  CORNERS  NEEDED  TO  GENERATE  THE 
DOMINATED  REGION  FOR  TIME  PERIOD  IJ  J<«ID(I) 
DIMENSION  I X < 1 5 » 25 » 2 ) » I D  ( 1 5  > »K<2) » J<3) »LL(25r2) 

NF=NT-1 

FIND  INTERSECTION  OF  DOMINATED  REGIONS 
DO  100  N=1 t NF 
Nl-N+1 
K(1)=ID(N) 

K<2)-ID(N1) 

J<1)=1 
J(2)=l 
J ( 3 ) =0 
L=1 

DETERMINE  MAXIMUM  VALUE  OF  C 

IF<IX<N»J<1)»1)*LT*IX<N1»J<2)»1)>  L=2 

M=N+L~1 

Ml-Nl+l-L 

SAVE  MAXIMUM  VALUE 
5  IXT=IX<M» J<L) » 1 ) 

L1-3-L 

10  J(L1)=J(L1)+1 

IF ( J ( L 1 ) . GT . K < L 1 ) )  GO  TO  15 
TEST  IF  DOMINATED  BY  NEXT  POINT 

IF(IX<M1*J<L1)»1) »LE. IXT  >  GO  TO  10 
15  J(L1)»J(L1)-1 
J  <  3 )  -  J  ( 3  >  H 
LL  <  J( 3) f 1 ) -IXT 

TF<IX(M»J<L)  » 2 )  ♦  LT  ♦  I X  <  M 1 »  J  v  L 1  >.?>>  L-L 1 
M-N+L-l 


L-L 1 


THIS  £0  VDO 

fSWk  COt  1  * 


1920 

M1=N1+1—L 

1930 

LL<  J<3>  >2>=IX<M»J(L)  *2) 

1940 

Ll=3-L 

1950 

IXT=IX(M1»J(L1) »2) 

1960 

20  J(L)=J(L)+1 

1970 

IF(J<L> .GT.K(L) >  GO  TO  30 

1980 

IF  ( I X  ( M  y  J  ( l_ )  y  2 ) .  L T  ♦  I XT  >  GO  TO  5 

1990 

J  ( 3 )  ~  J  ( 3 )  +  .1 

1995 

C 

NEW  CORNER 

2000 

LL  <J<3) » 1)  =  IX(M»J<L) »i) 

2010 

LL( J<3) »2)=IX<M» J(L) »2) 

2020 

GO  TO  20 

2030 

30  CONTINUE 

2040 

I D  ( N 1 )  -J  ( 3 ) 

2050 

DO  40  L~lrJ<3> 

2060 

IX<N1»L»1)=LL(L» 1 > 

2070 

40  I X ( N1 » L » 2 ) -LL ( L » 2  > 

2080 

100  CONTINUE 

2090 

RETURN 

2100 

END 

2110 

SIJBROUT I NE  DP  C  NT  >  IX*  I D  ?  SC ,  JS  ,  CS 1 1 P  > 

2111 

C 

THIS  SUBROUTINE  FINDS  THE  SOLUTION  TO  A  SUBPROBLEM 

2112 

C 

USING  DYNAMIC  PROGRAMMING 

2113 

C 

NT -THE  NUMBER  OF  TIME  PERIODS 

2114 

C 

IX ( I y  Jr  K) -SET  OF  POSSIBLE  SOLUTION  POINTS-  CORNERS 

2115 

c 

I -THE  TIME  PERIOD 

2116 

c 

J-THE  PAIR  INDEX 

2 .1 1 7 

c 

K-l  FOR  C  r 2  FOR  Y 

2118 

c 

I D <  I ) -  THE  SIZE  OF  IX  SET  FOR  TIME  PERIOD  1 1  ,.K-ID(  I ) 

2119 

0 

SC ( I r L ) "COSTS  IN  THE  OBJECTIVE  FUNCTION 

2119,1 

c 

I -THE  TIME  PERIOD 

2.119,2 

c 

L-l  TO  4 » ONLY  1  AMD  2  ARE  USED  IN  THE 

2119,3 

c 

OBJECTIVE  FUNCTION 

2119.4 

c 

J  S ( I y K ) = S  0  L U  T 1 0 N  P  OL ICY 

2119.5 

c 

I "TIME  PERIOD 

2119.6 

c 

K-l  FOR  C » 2  FOR  Y 

2119.7 

c 

JX ( I  •/ K  > -  THE  SOLUTION  POLICY  ASSUMING  A  PARTICULAR 

2119.8 

c 

POLICY  FOR  THIS  LAST  TIME  PERIOD  (I~NT> 

2119.9 

c 

I -THE  TIME  PERIOD 

2119.91 

c 

K-l  FOR  C  f 2  FOR  Y 

2120 

D I  MENS  I  ON  I X  ( 1 5 ,25r2>r  I D  <  1 5 ) » SC  ( 1 5  *  4  > » JS  (  i  '5,2),  JX  ( 1 5 » 2 ) 

2130 

CS-100000000. 

2140 

KF-ID ( NT ) 

2150 

IN- NT- 1 

2155 

c 

FIND  A  SOLUTION  FOR  EACH  POINT  IN  THE  LAST  TIME  PERIOD* NT > 

2160 

DO  1000  K-l f KF 

2170 

1 1-IX<  NT  f K , 1 > 

2180 

I2-TX ( NT r  K  y  2  > 

2185 

c 

DETERMINE  VALUE  OF  OBJECTIVE  FUNCTION  FOR  THE  LAST 

2186 

c 

TIME  PERIOD 

flUISPf.ei.  -  —  -  -  “*  ; 

(JRQW  06*  1  jfoivrtAbtuiv  iU  uX 
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2190 

2200 

2210 

2215 

2216 
2220 
2230 
2240 
2245 
2250 
2255 
2260 
2265 
2270 
2280 
2290 
2300 
2310 
2315 
2320 
2325 
2330 
2335 
2340 
2350 
2360 
2370 
2380 
2390 
2395 
2400 
2410 
2420 
2425 
2430 
2440 
2445 
2450 
2460 
2470 
2480 
2490 
2500 
2505 
250  , 
2510 
2520 
2530 
2540 
2550 


C=I1*SC(NT»1)+I2*SC(NT,2> 

JX  <  NT » 1 ) ”1 1 
JX  <  NT  *  2 )  =  12 

DETERMINE  BEST  POLICY  FOR  SUBSEQUENT  STAGES  (EARLIER 
TIME  PERIODS) 

DO  100  1  =  1  >  IN 
NI-NT-I 
L=0 

DETERMINE  OPTIMAL  Y 
10  L=L+1 

TEST  IF  <=  CONSTRAINT  IS  SATISFIED 
IF(IX(NI»L>2) »GT.I2>  GO  TO  10 
CALCULATE  CONTRIBUTION  TO  OBJECTIVE  FUNCTION 
CT= 1 2*SC ( N I  *  2 )  + 1 X ( N I >  L » 1 ) *SC < N I » 1 ) 

JC=IX<NI»L»1 > 

JY  =  I2 
L=L-1 
20  L=L+1 

DETERMINE  OPTIMAL  Y 

IF(L»GT . ID(NI ) )  GO  TO  50 
TEST  FOR  END  OF  STATES 

I F ( I X ( NI r  L  *  1 ) ♦ GT ♦ 1 1 )  GO  T050 
CALCULATE  CONTRIBUTION  TO  OBJECTIVE  FUNCTION 
CI“IX(NI »L»1>*SC<NI,1)+IX<NI,L»2>*SC<NI»2> 
IF(CI.GT.CT)  GO  TO  20 
CT=CI 

JC=IX(NIfL»l> 

JY=IX(NI »L*2) 

GO  TO  20 


TEST  FOR 


POLICY 


5 0  Cl "SC INI f 1 >*JX (NI  +  1 y 1 )+ S C ( N I » 2 > * J X < N I + 1 » 2  ) 

IF (Cl .GT.CT)  GO  TO  60 
CT  "Cl 

SAVE  POLICY 

JC~JX(NIH  r  1  ) 

JY"JX(NI  +  1  ?2> 

SAVE  OPTIMAL  STATES 
60  JX<NI»1)=JC 
JX  <  NI *  2 ) ~ JY 
C=C+CT 
Il-JC 
I2=JY 

100  CONTINUE 

TEST  EACH  SOLUTION  ASSOCIATED  WITH  A  STATE  OF  THE  LAST 
TIME  PERIOD  TO  DETERMINE  THE  MINIMUM  OBJECTIVE  FUNCTION 
IFIIP.EO.i  )  WRITE ( 6 r 901  )  (  C  t  JX  ( I  *  I.  )  r  JX<  I  *2)  »  J-t»  NT  > 
901  FORMAT  < FI 0.2, 21 10) 

IF(C.GT.CS)  GO  TO  1000 
CS=C 

DO  150  1=1 r NT 

tw  'i 

jtt  VUH 


t 

-  .  J| 

98  f 

2555 

C 

SAVE  SOLUTION 

2560 

JS ( 1 1 1 ) ~ JX (I , 1 ) 

2570 

JS ( I v  2 ) ~JX (1,2) 

2580 

150  CONTINUE 

2590 

1000  CONTINUE 

2600 

RETURN 

2610 

END 

2620 

SUBROUT  I NE  DUE  ( RHO ,  M » A  V » I C ,  I V ,  I R ,  IJ ,  REAR ) 

2621 

C 

THIS  SUBROUTINE  CALCULATES  THE  AVAILABILITY  CONSTRAINT 

2621.1 

C 

AND  R  BAR  (MEAN  NUMBER  OF  MACHINES  REPAIRED  PER  DAY) 

2621.2 

C 

IT  DETERMINES  THE  MINIMUM  VALUES  SATISFYING  THE  AVAILABILIT 

Y 

2621*3 

C 

CONSTRAINT ,  OR  IF  A  (C,Y)  PAIR  IS  SPECIFIED ,  DETERMINES 

2621.4 

C 

IE  IT  SATISFIES  THE  CONSTRAINT. 

2621.5 

C 

RHQ-RHO  (LAMBDA  BAR  X  M/MU) 

2621.6 

C 

M~THE  NUMBER  OF  MACHINES 

2621 .7 

C 

AV-THE  AVAILABILITY  REQUIREMENT 

2621.8 

C 

IOC 

2621.9 

C 

I  Y~Y 

2622 

C 

IR=S«ITCHi»i  TEST  IF  AVAILABILITY  CONSTRAINT  IS 

2622.1 

C 

SATISFIED  FOR  A  SPECIFIED  (C,Y>  PAIR 

2622 . 2 

C 

NOT  1  GENERATE  (CvY)  PAIR 

2622.3 

C 

U~MU  (REPAIR  RATE) 

2622 . 4 

c 

RBAR-RBAR <  MEAN  NUMBER  OF  MACHINES  REPAIRED) 

2625 

C**###CALCULATE  A  1 

2630 

A  AV/  ( 1  ~AV> 

2635 

C 

TEST  IF  OY 

2640 

IF(IC.GT.IY)  GOTO  600 

2645 

C 

TEST  IF  O  RHO 

2650 

10  RORHO/IC 

< 

2660 

I F  ( RC .  NE .  1 )  RR-- 1  /  ( .1.  -  RC ) 

f 

r 

2670 

S2:"0 

2680 

T  2  ■■■■•■  :l. 

2685 

C 

CALCULATE  S2 

2690 

DO  100  J-"  1 ,  M 

2700 

82~S2+T2 

2710 

IFCT2.LT. .00000001 )  T2«0 

2720 

100  T2-™T2#RH0/(M*IC)#(M-J) 

2730 

IF(IR.EQ.l)  GO  TO  150 

2740 

IF(RC.LE.l)  GOTO  150 

2745 

C 

TEST  IF  AVAILABILITY  CONSTRAIN  CAN  BE  SATISFIED  FOR  Y>C 

2750 

IF(A*S2.LE.~RR>  GOTO  150 

2760 

IC---TC+J 

2770 

GO  TO  10 

2780 

150  Sl~(> 

2785 

c: 

CALCULATE  31 

2790 

reel. 

2800 

DO  200  T  ~  1  *  fC 

2810 

Sl=si+Tl 

2820 

I F  (T1.IJ..  0 0  0 0 0 0 0 .1  )  r  1  •••*)  j 

'  i 

i 

ENOW  COil  1  -  •  •  i 

L. 

_ 1 

99 


i 


2825 

C 

CALCULATE  Y 

2830 

200 

T1=T1*RH0/I 

2840 

IF C IR.EQ. 1 )  GOTO  400 

2850 

TF  =A*S2-S1/T1 

2860 

IFCRC.NE.l)  TF-ALOGC <S1/T1+PR>/CA*S2+RR>  >/ALOGCRC 

2870 

IY=IC+TF+. 999999 

2880 

IFCTF.LT*-. 0001)  GO  TO  *^00 

2890 

GO  TO  1000  m 

2900 

400 

T3=RC**CIY-IC> 

2905 

C 

TEST  IF  CC*Y)  SATISFY  AVAILABILITY  CONSTRAINT 

2910 

S1=S1+T1*C1.~T3)*RR 

2920 

T1=T1*T3  1 

2930 

S2=T1*S2 

2940 

GO  TO  2000 

2945 

C 

Y<C 

% 

2950 

500 

IY=IC 

2960 

600 

S1=0. 

2970 

T 1  - 1 . 

2975 

C 

CALCULATE  AVAILABILITY  CONSTRAINT 

2980 

DO  700  1*1, IY 

2990 

S1-S1+T1 

3000 

IFCT1.LT.. 00000001)  T1=0 

3010 

700 

T1=T1#RH0/I 

3020 

S2~0 . 

3030 

T2*l. 

3040 

DO  800  J~1 » IC--I Y 

3050 

S2=S2+T2 

3060 

IFCT2.LT. .00000001)  T2-© 

3070 

800 

T2=T2/M*RH0/ C I Y+ J  >  *  C  M-J ) 

3080 

DO  900  J=IC-IY+1»M* 

3090 

S2=S2+T2 

3100 

IFCT2.LT. .00000001 )  T2=0 

3110 

900 

T2=T2/IC*RH0/M*CM-J> 

3120 

S2=S2#T1 

3125 

C 

ALREADY  KNOW  CC,Y).  WANT  TO  TEST  IF  AVAILABILITY 

3126 

C 

CONSTRAINT  IS  SATISFIED  AND  WANT  TO  CALCULATE  RBAR 

3130 

IFC IR.EQ. 1)  GO  TO  2000 

3135 

C 

TEST  IF  AVAILABILITY  CONSTRAINT  IS  SATISFIED 

3140 

IF  CS1.GE.A*S2>  GO  TO  1000 

3145 

C 

INCREASE  C 

3150 

IOXC+1 

3155 

C 

TEST  IF  OY+M 

3160 

TFC  IC.LE.IY+M)  GO  TO  600 

3170 

1000 

RETURN 

3175 

C 

CALCULATE  RBAR 

3180 

2000 

RBAR-RHO/U*  C  Si  +S2  )  /  C  SI  +RH0/.IC*S2+T1 ) 

3185 

C 

AVAILABILITY  CONSTRAINT  NOT  SATISFIED 

3190 

IF  CS1.LT.A*S2>  IR--1  t 

3200 

RETURN  • 

3210 

END 

1 
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APPENDIX  B 


* 


ALGEBRAIC  REFORMULATIONS  OF 
THE  AVAILABILITY  CONSTRAINT 


This  appendix  will  give  the  reconfiguration  of  the  avail¬ 
ability  constraint  in  a  more  convenient  form.  The  availability 
constraint  is  given  by: 

q  >  0.90  . 

n  == 

This  can  be  generalized  for  any  availability  to 

qn  =  a  * 

where  a  is  the  availability. 


£ 

n=0 


Using  the  expression  (II-4)  for  q  (and  noting  that 

n 


n  <  y-1)  ,  we  have 


D  Mpr 


n=0 


y+M 

M  -  £  (W)p £ 

l=y 


>  a  . 


This  can  be  rewritten  as 


£  pr 

n=0 


1  '  j?o  «  p3+y 


>  a  . 


100 


101 


Finally, 


% p"  4  (*)  %  fcK* 


(B-l) 
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